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ABSTRACT 

To  provide  an  understanding  of  the  acoustic  source  characteristics  of  the 
boundary  layer  transition  r^on,  unsteady  velocity  field  measuremrats  of  an  isolated, 
artificially  generated  turbulent  spot  were  made  in  a  zmt>-pressuFe  gradient  laminar 
boundary  layer.  These  measurements  were  performed  in  a  water  channel,  using  a 
laser-Doppler  velocimeter.  They  provide  quantitative  information  describing  the 
large-scale  unsteady  displacement  thickness  fluctuations  due  to  the  passage  of  a 
turbulent  spot.  Fluctuations  of  the  displacemrat  thickness  are  related  to  the  radiated 
noise  through  the  Liq)mann  acoustic  analogy.  The  description  of  the  large-scale 
velocity  field  is  based  on  ensemble-averaged  unsteady  velocity  data,  from  which  the 
displacement  thickness  is  calculated.  These  results  are  used  to  calculate  the  velocity 
normal  to  the  plate,  as  well  as  the  characteristic  rise  times,  tj,  of  the  di^lacement 
thickness.  Comparison  of  the  results  at  diffnent  streamwise  stations  shows  that  the 
rise  time  and  the  mass  flux  deficit  peak  amplitudes  increase  with  distance  from  the 
generation  point,  while  the  peak  amplitudes  of  the  normal  velocity  increase  early  in 
the  spot  development,  but  level  off.  For  all  these  parameters,  the  rate  of  spatial 
growth  is  greatest  in  the  stations  closest  to  the  generation  point.  The  normalized  rise 
time,  tjUe/Ax  was  found  to  span  a  range  of  values  0.06  <tiUe/Ax<  0.74  over  the 
locations  measured,  where  Ue  is  the  spot  convection  speed  and  Ax  is  a  transition  zone 
Iragth  obtained  fiom  a  correlation.  The  lowest  values  of  this  paramet^  were  found 
to  exist  in  the  upstream  end  of  the  transition  zone.  A  scaling  analysis  using  these 
experimental  results  suggests  that  the  sound  radiated  by  the  large-scale  motion  due  to 
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an  isolated  turbulent  spot  increases  as  the  spot  grows  and  has  a  dipole  character. 
Extrapolation  to  a  natural  transition  zone  indicates  that  sound  radiation  from  the  large- 
scale  intermittent  motion  is  highest  in  the  middle  part  of  the  transition  zone  because 
the  spot  density  is  highest  there.  To  estimate  the  effect  of  turbulent  spot  interaction 
on  the  sound  radiation  from  a  single  spot,  measurements  were  carried  out  for  two 
turbulent  spots  generated  at  the  same  location  with  a  fixed  time  delay.  The  second 
spot  partially  merged  with  the  first  at  the  furthest  downstream  station.  The  spot 
interaction  was  not  found  to  alter  the  isolated  spot  results  significantly. 
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CHAPTER  1 


INTRODUCTION 


1.1  OVERVIEW 

Although  the  production  of  sound  by  unsteady  fluid  modem  has  been  a  topic 
of  interest  for  nearly  four  decades,  the  pioneering  work  of  Lighthill  (1952,  1954) 
remains  the  point  of  reference  for  much  of  our  currrat  understanding  of  the  subject. 
The  focus  of  Lighthill’s  early  work  in  this  area  was  the  production  of  sound  by  jets. 
He  was  motivated  by  the  needs  of  jet  propulsion  which  was  in  its  early  stages.  It  was 
realized  that  the  presence  of  solid  bodies  could  substantially  alter,  even  increase,  the 
aerodynamic  productiem  of  sound  (see,  e.g. ,  Curie,  1955).  Sound  generated  by  a  flow 
interacting  with  a  solid  surface  may  be  classified  into  two  cat^ories  according  to  tu>w 
the  flow  unsteadiness  is  generated.  If  the  dominant  noise  generation  mechanism  is 
unsteady  flow  impinging  on  a  solid  body,  the  sound  produced  is  refinred  to  as 
"interaction  noise."  This  type  of  noise  production  occurs  when  helicopter  blades 
interact  with  the  blade-tip  vortices  th^  dmy  generate  and  in  a  jet  engine  compressor 
stage  due  to  rotor/stator  interaction.  When  the  unsteady  flow  responsible  for  the 
sound  productiim  results  from  viscous  action  at  a  solid  body  it  is  called  "self-noise." 
Examples  include  sound  ftom  a  turbulent  boundary  layer  on  a  wall  (v  from  the  wake 
of  a  circular  cylinder. 

Another  mcample  of  "self-noise”  is  the  noise  produced  by  a  boundary  layo* 
undergoing  transition  from  laminar  to  turbulent  flow.  The  principle  noise  generation 
is  thought  to  occur  in  the  "intermittoit  r^km,”  where  the  final  stage  of  the  transition 
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between  laminar  and  turbulent  flow  causes  large  fluctuations  in  the  velocity  and 
pressure  fields.  Not  only  are  these  fluctuations  large,  but  they  are  thought  to  occur 
such  that  the  flow  is  an  efficient  acoustic  radiator  (Lauchle,  1981,  1989;  Somette  and 
Lagier,  1984a,  b;  Lagier  and  Somette,  1986;  and  Audet  et  al.  1989a,  b).  Lauchle 
(1981),  Somette  and  Lagier  (1984a,  b),  Lagier  and  Somette  (1986),  and  Audet  et 
a/.  (1989)  predicted  a  monopole  type  source.  Lauchle  (1989)  showed  later,  using 
scaling  arguments,  that  a  naturally  transitioning  boundary  layer  radiates  sound  power 
proportional  to  only  slightly  more  efficiently  than  a  fully  turbulent  boundary 
layer.  In  flows  of  low  Mach  number,  a  monopole  source  is  the  most  efficient  at 
radiating  sound  and  would  therefore  be  expected  to  dominate  the  generated  sound 
field.  If  boundary  layer  transition  is  an  efficient  sound  source,  the  importance  to  the 
silencing  of  underwater  vehicles  is  immediately  apparent. 

There  is  some  experimental  evidence  which  suggests  the  predominance  of 
transition  as  a  noise  source.  Perraud  (1989)  conducted  measurements  of  wall  pressure 
on  a  body  of  revolution  in  a  wind  tunnel.  By  placing  wall-mounted  microphones  in 
the  laminar  portion  of  the  boundary  layer  upstream  of  the  transition  region  and  by 
using  conditional  sampling  triggered  on  ^t  arrivals  at  a  downstream  location  he 
showed  that  the  spots  did  radiate  sound  at  ^xpreciable  levels.  Perraud  further  showed 
that  the  limits  set  by  Lauchle  (1981)  for  the  sound  spectral  levels  bounded  his  own 
wall-pressure  spectra.  Attempts  to  rq)eat  these  experiments  in  a  water  towing  tank 
were  unsuccessful  because  of  the  mechanical  noise  from  the  towing  mechanism. 
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Audet  et  a/.  (1989a,  b)  also  demonstrated  that  transitional  flow  is  an  ai^reciable 
noise  source.  This  work  was  a  combined  analytical/experimental  study.  Using  a 
model  based  on  the  previous  work  of  Somette  and  Lagier  (1984a,  b)  and  Lagier  and 
Somette  (1986),  they  clearly  sqiarated  the  waU  pressure  fluctuations  into  contributions 
due  to  (1)  the  effect  of  the  turbulence  in  the  turbulent  spots  and  (2)  the  radiated  sound 
caused  by  the  large-scale  fluctuations  in  the  boundary  layer  flowfield.  The  analysis 
showed  an  s^reciable  contribution  to  the  wall-pressure  due  to  pressure  fluctuations 
in  the  acoustic  wavoiumber  domain.  The  confirmation  of  these  results  consisted  of 
wall-pressure  measurements  in  a  transitional  boundary  layer  which  were  corrected  for 
the  spatial  filtering  effects  of  the  microiAone  size.  These  corrected  RMS  pressure 
measurements  were  shown  to  give  good  agreement  with  the  sum  of  the  contributions 
of  the  two  parts  of  the  wall-pressure,  the  acoustic  contribution  due  to  the 
intermittency,  and  the  intermittency-weighted  turbulent  source. 


1.2  BOUNDARY  LAYER  TRANSITION  FLOWS 

A  full  description  of  the  noise  from  a  naturally  occurring  boundary  layer 
transition  requires  a  knowledge  of  the  transition  flowfield.  Since  this  process  is  to 
some  extent  a  "random"  one,  its  descripticm  is  generaUy  given  in  statistical  trains.  A 
laminar  boundary  layer  may  undergo  a  transition  to  turbulent  flow  by  many  routes. 
The  classical  scenario  begins  with  the  laminar  boundary  layer  exhibiting  two- 
dimraisional  viscous  ToUmirai-Schlichting  instability  waves  that  grow  exponentially 
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from  their  inception.  This  instability  mode  itself  becomes  unstable  and  is  susceptible 
to  disturbances  with  spanwise  nonuniformity  and  thus  it  begins  to  develop  a  spanwise 
waviness.  Up  to  this  point  the  energy  contained  in  the  instability  wave  fluctuations 
is  only  a  very  small  lotion  of  the  total  flow  oiergy.  The  three-dimensional 
instability  waves  thra  interact  to  form  local  high-shear  regions  which  cause  the 
laminar  flow  to  break  down  rapidly  into  turbulence.  As  these  localized  regions  of 
turbulence,  or  turbulent  spots,  convect  downstream,  they  grow  by  destabilizing  the 
surrounding  laminar  boundary  layer  (Gad-el-Hak,  et  al. ,  1981),  causing  it  to  become 
turbulent  at  the  edges  of  the  spot.  The  spots  grow  until  they  merge  with  each  other 
and  a  turbulent  boundary  layer  is  formed.  The  streamwise  extent  of  this  region  is 
small  compared  to  the  other  regions  of  the  flow. 

The  transition  process  need  not  proceed  by  this  path,  however.  The  process 
described  above  is  based  on  the  behavior  of  an  unstable  laminar  boundary  layer, 
driven  by  a  low-turbulence  outer  flow,  in  the  presence  of  infinitesimal  disturbances. 
If  the  freestream  turbulence  level  is  high,  or  if  there  is  some  large  disturbance  to  the 
boundary  layer,  then  the  early  stages  of  transition  may  be  bypassed.  An  example  is 
a  laminar  boundary  layer  separation  in  which  the  laminar  shear  layer  becomes 
turbulent  before  re-attaching  to  the  body.  A  boundary  layer  subjected  to  an  extreme 
adverse  pressure  gradient  may  also  bypass  the  earlier  stages  of  transition,  even 
without  sqxiration. 

By  whatever  means  the  final  stage  of  transition  is  reached,  it  is  characterized 
by  its  intermittent  nature.  Early  studies  of  the  transition  zone  (Dhawan  and 
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Narasimha,  19S8)  characterized  it  as  a  region  alternating  locally  between  a  fiilly- 
turbulent  and  a  fiilly-laminar  structure.  By  defining  an  intermittency  parameter,  7, 
as  the  percentage  of  the  time  that  the  flow  at  a  given  point  is  turbulent,  the  time- 
averaged  velocity  profile  and  its  int^ral  parameters,  such  as  di^lacement  thickness, 
may  be  described  using  an  intermittency-weighted  form; 

V(%)  "  (1  -  7(i))5\(x)  ^  y(x)8-(x)„  (1-» 

where  S\  is  the  fiilly-laminar  (Kasius)  displacement  thickness  and  6\  is  the  fully 
turbulent  time-averaged  displacement  thickness.  The  intermittency  parameter  has  beoi 
shown  to  represent  adequately  the  avmge  properties  of  the  flowfleld  for  a  zero- 
pressure  gradient  boundary  layer  transition.  For  flows  with  a  significant  pressure 
gradirat,  there  does  not  appear  to  be  such  a  universal  scaling  parameter  (Narasimha, 
etaL,  1984). 

The  intermittent  nature  of  the  flow  is  shown  in  the  ensemble-averaged  velocity 
profile  measurements  of  Wygnanski  et  a/. (1976).  The  averaged  profile  is  laminar 
upstream  of  the  spot,  but  as  the  spot  passes  by,  the  near-wall  regions  of  the  flow 
undergo  a  rapid  acceleration,  and  then  a  deceleration,  corre^nding  to  a  change  from 
a  laminar  velocity  profile  to  a  turbuleit  i^file  and  back  again.  This  eveit  is 
accompanied  by  an  injection  of  high-momeitum  fluid  from  the  fleestream,  followed 
by  an  ejection  of  the  low-momentum  fluid  flom  near  the  wall  as  the  flowfleld  returns 
to  its  laminar  state.  The  final  return  to  the  initial  laminar  state  proceeds  very  slowly 
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and  is  characterized  by  a  laminar  "calming  region”  (Schubauer  and  Klebanoff,  1955). 

The  alternating  appearance  of  laminar  and  turbulent  flow  causes  velocity  fluctuations 
normal  to  the  plate  which  are  larger  in  magnitude  than  those  occurring  in  a  fully 
turbulent  layer. 

The  large-scale  motion  exhibited  by  an  ensemble-averaged  rq>resentation  of 
the  spot  initially  led  some  to  view  the  i^t  as  a  large,  arrowhead-shaped  vortex  tube 
(Coles  and  Barker,  1975;  Cantwell  et  al.,  1978).  This  was  shown  to  be  an 
oversimplification  by  the  flow  visualization  studies  of  Cantwell  et  al.  (1978);  Matsui 
(1980);  Gad-el-Hak  et  al.  (1981);  and  Perry  et  al.  (1981),  which  indicated  that  no 
such  large-scale  structures  existed  ~  instead,  the  spot  is  composed  of  small  structures 
not  unlike  those  found  in  a  fiilly  turbuloit  boundary  layer.  Later  investigations  (see, 
e.g.,  Wygnanski  et  al.,  1981;  Sankaran  et  al.,  1991)  demonstrated  that  the  ensemble 
averaged  flowfield  is  a  somewhat  misleading  rq>resratation  of  the  unsteady  flow  in 
a  turbulent  spot.  However,  Van  Atta  et  a/.  (1982)  showed  that  the  potential  flow 
velocity  responds  to  a  ^t  passage  as  it  would  to  a  "bump”  in  a  wall  described  by 
y=d*(x,t),  where  5*  is  the  local  displsu^ment  thickness  found  from  the  ensemble- 
averaged  unsteady  velocity  field.  Furthermore,  the  scaling  results  obtained  from  an 
ensemble-averaged  velocity  field  do  provide  accurate  measurements  of  the  ovoall 
scales  of  the  spot  such  as  the  streamwise,  spanwise,  and  normal  growth  rates,  as  well 
as  the  convection  velocity  of  the  spot  (Sankaran,  et  al.,  1991).  Thus,  although  the 
rasemble-averaged  description  of  the  velocity  fleld  is  somewhat  misleading  in 
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describing  the  details  of  the  boundary  layer  flow,  it  may  be  used  to  predict  the  effect 
of  the  unsteady  boundary  layer  on  the  potential  flow. 

The  displacement  effect  is  an  important  one  in  considering  sound  goierated  by 
a  viscous  flow  (Liepmann,  1954).  It  is  wdl  known  that  the  effect  of  a  boundary  layo* 
on  the  surrounding  potential  flow  is  equivalent  to  an  inviscid  flow  over  a  wall  of 
height  equal  to  6*.  The  displacement  thickness  is  an  integral  paramet^  expressing  the 
amount  of  mass  flux  lost  to  an  incoming  stream  due  to  the  presence  of  the  boundary 
layer.  The  mass  flux  lost  to  the  incoming  stream  may  be  directed  either  in  the 
spanwise  direction  or  away  from  the  solid  surface.  Another  important 
consideration  is  the  range  of  temporal  and  spatial  scales  of  the  unsteady  flow.  In  a 
naturally-occurring  transition  there  are  three  length  scales.  The  first  is  the 
charactmistic  size  of  the  turbuloit  eddies  of  which  the  turbulent  spots  are  composed. 
The  next  level  is  the  size  of  the  spots  themselves.  Lastly,  there  is  the  streamwise 
extent  of  the  entire  transition  zone.  Similarly,  temporal  scales  exist  at  these  same 
levels:  the  time  for  a  turbuloit  eddy  to  pass  a  given  location,  the  passage  time  for  a 
turbulent  spot,  and  the  time  between  spot  passages.  The  transition  zone  involves  a 
range  of  scales  that  are  larger  than  those  found  in  a  fully  turbulrat  boundary  layer. 
These  scales  can  provide  oider-of-magnitude  estimates  for  the  frequencies  of  the 
sound  that  will  be  present  in  the  furfield. 
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1.3.  ACOUSTICS  OF  BOUNDARY  LAYER  TRANSITION 

Given  that  any  unsteady  flow  has  the  potential  to  cause  sound  goieration,  it  is 
apparent  that  the  various  transitional  flow  structures  are  possible  sources.  The 
important  distinction  between  the  tranational  and  fully  turbulent  portions  of  the  flow 
is  that  the  transitional  regimes  are  characterized  by  boundaries  between  the  laminar 
(steady)  and  turbulent  (unsteady)  regions  of  the  boundary  layer.  Because  the  volume 
of  the  unsteady  regions  grows  as  they  convect  downstream,  the  acoustic  "source 
volume"  increases  with  time.  In  the  intermittoit  region,  then,  there  is  not  only  the 
mechanism  by  which  a  fiiUy  turbuloit  layo*  would  radiate  sound,  but  also  an 
additional  effect  due  to  the  localization  and  growth  of  the  turbulent  regions  (Lauchle, 
1980,  1981;  Lagier  and  Somette,  1984a,  1984b,  1986). 

The  problem  of  noise  from  a  transitioning  boundary  layer  did  not  receive  much 
attention  until  the  early  1980s.  Most  of  die  work  on  boundary  layer  noise  has  dealt 
with  radiation  from  fully  turbulent  boundary  layers.  Sound  radiation  from  a  boundary 
layer  was  studied  by  Phillips  (1955)  and  Powell  (1960).  Most  of  this  early  work  used 
Curie’s  (1955)  extrasion  of  Lighthill’s  (1952,  1954)  analysis  to  include  the  effects  of 
solid  bodies.  Specifically,  it  dealt  with  how  to  describe  the  acoustic  sources  in  toms 
of  distributions  both  of  quadrupoles  due  to  unsteady  normal  stresses  in  the  boundary 
layer  and  of  dipoles  due  to  unsteady  forces  at  the  solid  wall. 

For  a  large  plate  ("large"  being  defined  as  the  length  scale  of  the  plate  bdng 
much  larger  than  the  acoustic  waveloigdi  of  intoest)  Powell  (1960)  showed  that  the 
effect  of  the  solid  wall  was  to  reflect  the  sources  in  the  flow.  He  treated  this 
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reflection  by  replacing  the  wall  with  a  set  of  sources  of  equal  strength  but  positioned 
so  that  there  is  a  symmetry  between  a  source  and  its  opposite  virtual  number.  The 
net  effect  of  this  reflection  is  that  lateral  quadrupoles  become  octupoles,  longitudinal 
quadrupoles  double  in  strength,  and  surface  dipoles  become  quadrupoles.  Based  on 
this  analysis  the  radiation  from  a  turbul«it  boundary  layer  should  have  at  most 
quadrupolar  radiation  characteristics. 

The  classification  of  aeroacoustic  sources  into  monopoles,  dipoles,  etc.,  is 
useful  if  handled  with  care  (Powell,  1990).  If  such  a  division  of  sources  is 
meaningful,  it  can  provide  a  straightforward  means  of  ranking  source  mechanisms  in 
toms  of  their  efficiency,  especially  for  flows  at  low  Mach  number.  The  effidoicy 
of  a  compact  monopole  scales  with  Mach  number,  a  dipole  with  Mach  number  cubed, 
a  quadrupole  with  Mach  number  to  the  fifth  power.  For  underwater  flows  this  is 
especially  useful,  since  the  Mach  number  seldom  goes  above  0.01 .  This  explains  why 
so  much  effort  has  been  expended  in  determining  how  the  dipoles  due  to  wall  pressure 
fluctuations  manifest  themselves  in  turbulent  boundary  layer  noise  [see,  e.g.,  Phillips 
(19SS);  Tam  (1975);  and  Hardin  (1991)].  If  wall  pressure  fluctuations  bdiave  as 
dipole  sources,  they  would  dominate  the  radiation  field  of  a  low  Mach  number  flow 
and  attention  can  be  focussed  on  that  particular  mechanism. 

Tam  (1975)  used  a  differait  sq^proach  to  the  study  of  turbulent  boundary  layer 
noise.  Using  an  empirical  wall-pressure  q)ectrum  developed  by  Maestrello  (1965, 
1967)  he  estimated  the  pressure  fluctuations  on  a  plane  parallel  to  the  flat  plate,  but 
just  above  the  boundary  layer.  The  pressure  on  this  surface  was  then  related  to  the 
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pressure  at  a  point  in  the  farfield  using  the  method  of  stationary  phase.  This  approach 
allowed  Tam  to  show  that  the  portion  of  the  source  spectrum  that  radiates  sound  is 
that  for  which  the  wavenumber  components  have  supersonic  phase  speeds.  Continuing 
this  iq)proach,  Hardin  (1991)  demonstrated  that  localized  structures  are  more  likely 
to  produce  this  kind  of  spectrum,  precisely  because  of  their  localization  in  space. 
From  this  type  of  argument  Hardin  wait  on  to  hypothesize  that  in  a  turbulent 
boundary  layer  horseshoe  vortices  are  the  primary  source  of  sound  radiadcm. 

Liepmann  (1954)  [see  also  Laufo*  et  n/.(1964)]  had  already  anticipated  such 
an  s^roach.  His  particular  method  consisted  of  relating  the  radiation  from  a 
boundary  layer  to  fluctuations  in  the  velocity  normal  to  the  solid  surfsure,  which  may 
be  deduced  by  measuring  the  boundary  layer  displacement  thickness,  5*. 

If  the  boundary  layer  is  unsteady,  then  it  is  expected  that  the  mass  flux  normal 
to  the  wall  is  also  unsteady  and  that  these  fluctuations  give  rise  to  sound  radiation  if 
they  contain  frequency  components  that  possess  supersonic  phase  qieeds.  The 
acoustic  problem  is  then  solved  by  specifying  the  behavior  of  the  di^lacement 
thickness  fluctuations  from  the  known  boundary  layer  velocity  field  and  using  the 
deduced  normal  velocity  distribution  as  a  boundary  condition  for  the  acoustic  wave 
equation. 

This  situation  is  similar  to  the  flow  over  a  wavy  wall  of  hdght  y=d*,  as 
described  by  Van  Atta  et  a/.  (1982).  The  radiation  from  a  compact  flow  ova*  a 
compact  surface  is  given  by  Crighton  (1975).  The  dominant  source  is  a  nranopole 
caused  by  the  rate  of  change  at  which  the  volume  of  fluid  is  diq>laced  by  the  wall. 
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Crighton  (1975)  points  out,  however,  that  this  source  mechanism  exists  only  if  the 
there  is  a  net  unsteady  mass  flux  over  the  oitire  surface.  If  there  is  mass  flux  away 
from  the  wall  in  one  location,  but  an  equal  amount  towards  the  plate  in  another 
location,  there  is  no  net  mass  flux  and  there  is  no  monopole  radiation  field  as  Icmg  as 
the  positive  and  negative  peaks  are  in  close  enough  proximity  to  be  considered  a 
compact  source  distribution. 

A  physical  system  which  behaves  like  a  monopole  is  a  ^here  that  undergoes 
spherically  symmetric  oscillations  of  its  size,  i.e.,  the  motion  of  its  surface  is  tmly 
radial  (see  Figure  1 .  la).  The  q)here  di^laces  the  fluid  uniformly  in  an  amount  equal 
to  the  rate  of  change  of  its  volume.  In  addition,  there  is  a  net  mass  flux  in  the  radial 
direction  at  any  given  time  when  the  ^)hoe  is  eithor  expanding  or  contracting.  For 
higher  order  sources  the  instantaneous  net  mass  flux  normal  to  the  surface  is 
identically  zero.  For  example,  a  sphme  whose  surface  deforms  in  a  non-axially 
symmetric  manner,  but  whose  volume  is  ccmsmved  (see  Figure  1.1b),  radiates  like  a 
quadrupole.  While  the  qntial  distribution  of  the  qrfiere’s  volume  changes  with  time, 
the  total  volume  does  not.  These  non-symmetric  spatial  variations  cause  {riiase 
cancellations  of  the  pressure  which  lead  to  a  lower  radiation  efficiency.  A  boundary 
layer  will  produce  monopole  radiation  if  the  volume,  whose  bounding  surfaces  are  the 
solid  wall  and  die  sut^kx  y  -  £*(x,t),  changes  with  time. 

The  Uqimann  sqiproach  has  not  recdved  as  much  attention  as  those  based  on 
a  Lighthill  formulation.  Lauchle  (1981),  however,  used  the  Liqimann  analogy  to 
describe  the  radiation  from  the  intomittent  zone  of  a  flat  plate  boundary  layer.  He 
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^  [a(t)  •  n]d§  =  g(t)  ^  0 
S 


(a)  Small  radially  oscillating  sphere  models  a  numopole.  A  mmiopole  is 
distinguished  by  a  net  volume  displacement,  reflected  in  a  nonz^  unsteady  mass 
flux  normal  to  the  surface  S  surrounding  the  ^hete. 


at 

^[ur(t) -nJdS  =0 


(b)  A  quadrupole  may  be  moddled  by  a  ^>here  undogoing  elliptical  distortions. 
The  quadrupole  is  distinguished  by  its  consmvation  of  volume,  reflected  in  die 
constant  zero  net  mass  flux  in  the  radial  direction. 


Figure  1.1  Basic  source  types  and  thdr  normal  mass  flux  distributions. 
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hypothesized  that  the  rapid  fluctuations  in  displacement  thickness  at  the  edges  of 
turbulent  spots,  which  could  be  characterized  by  a  time  scale  tj,  were  a  weak 
monopole  source  mechanism.  Lacking  information  regarding  the  behavior  of  these 
fluctuations,  he  wrote  an  expression  for  the  assumed  behavior  using  a  conditioned- 
signal  approach  (Libby,  1975)  and  assumed  an  intermittency-weighted  form  for  the 
unsteady  displacement  thickness.  Lauchle  (1981)  estimated  the  radiation  intensity 
(power  per  unit  area)  and  found  it  to  be  significantly  greater  than  that  radiated  from 
a  fully  turbulent  boundary  layer  covering  the  same  surface  area.  He  also  found  that 
the  high  frequency  portions  and  peak  values  of  the  calculated  farfield  spectra  were 
very  sensitive  to  the  value  of  tj. 

Later  work  by  Lauchle  (1989)  modified  his  earlier  conclusions  and  indicated 
that  the  intermittent  flow  is  a  monopole  type  source  only  for  low  frequencies. 
Another  way  to  state  this  is  that  the  source  is  no  longer  compact  at  high  frequencies 
and  as  a  result  a  classical  aeroacousdcs  treatment  does  not  necessarily  apply.  As 
stated  above,  a  boundary  layer  producing  monopole-like  radiation  will  exhibit  the 
following  bdiavior:  (1)  the  volume  bounded  by  y=5*  and  y=0  should  change  with 
time  and  (2)  the  unsteady  net  mass  flux  normal  to  the  wall  is  not  identically  zero. 
These  requirements  imply  either  that  turbulent  spots  must  grow  in  size  or  that  the 
number  of  spots  must  increase  in  the  downstream  direction  in  order  to  radiate  noise. 
Therefore,  the  monopole  nature  of  the  source  may  be  determined  by  specifying  the 
behavior  of  the  unsteady  mass  flux  normal  to  the  wall  as  a  function  of  position  in  the 
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intermittent  region.  If  the  normal  mass  flux  is  both  unsteady  and  nonzero  then  a 
monopole  does  exist  for  low  frequencies. 

Another  tq^proach  to  the  analysis  of  transition  noise  was  that  of  Lagier  and 
Somette  (1984a,  1984b,  1986).  They  hypothesized  that  in  intermittoit  flows  such  as 
boundary  layer  transition,  the  transioits  caused  by  the  intermittency  would  lead  to  an 
efficient  source  generation  mechanism.  Hiis  analysis  proceeded  by  making  use  of  an 
intomittency-weighted  form  of  the  lighthiU  acoustic  analogy.  The  intermittency- 
weighting  assumption  was  justified  by  considering  the  transition  region  to  be  a  mixture 
of  two  fluids,  one  laminar  and  one  turbulent.  The  form  of  the  acoustic  source  for 
boundary  layer  transition  determined  in  this  manner  was  found  to  include  the  large- 
scale  noise  source  modeled  by  Lauchle  (1981).  This  leading  term  was  interpreted  to 
be  a  monopole  caused  by  the  growth  of  the  turbulent  "phase”  of  the  flow  at  the 
expoise  of  the  laminar  "phase."  Hie  equivalrace  between  the  results  of  the  two 
sq/proaches  is  shown  explicitly  by  Audet  et  a/.  (1989a,  1989b). 

Previous  analyses  of  the  acoustic  generation  mechanisms  have  thus  relied  cm 
assumptions  concerning  the  behavior  of  the  boundary  layer  flowfield.  The  presort 
work  is  intended  to  supply  measurements  of  the  unsteady  velocity  tield  present  in  the 
principle  feature  of  transition,  the  turbulent  spot.  From  this  data  the  time  scale  1^  and 
the  acoustic  source  strength  for  a  turbulent  spot  is  estimated.  A  discussion  of  the 
implications  of  the  measured  flowfield  behavior  for  existing  models  of  transition 
noise,  e.g. ,  the  ranking  of  the  source  behavior  in  terms  of  monopoles,  dipoles,  etc. , 
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is  also  presented.  Finally,  the  data  is  used  to  provide  estimates  of  the  sound  radiated 
from  a  turbulent  spot  and  to  evaluate  the  assumptions  of  previous  investigators. 

1.4  PROJECT  GOALS 

The  fluid  mechanics  issues  that  must  be  addressed  in  order  to  further  the 
understanding  of  transition  noise  are: 

1.  The  proper  form  for  the  displacement  thickness  fluctuation  and  its 
relation  to  the  normal  velocity  derivative  fluctuation. 

2.  The  value  of  the  time  scale,  tj,  and  its  variation  with  position. 

3.  The  growth  rate  of  the  volume  of  sound-producing  fluid  (i.e.,  of  the 
turbulent  spot)  and  its  vnriatitm  with  position. 

4.  The  degree  of  attenuation  of  these  parameters  due  to  interaction 
between  turbuloit  spots. 

These  questions  have  an  important  bearing  on  the  relevant  acoustics  issues: 

1.  The  radiation  characteristics  of  a  turbulent  spot  at  different  stages  of  its 
development. 

2.  The  effect  of  the  form  of  the  displacement  thickness  fluctuation  on  the 
sound  radiation  estimate.  Here  the  int^est  is  in  what  effect  the 
simplifying  assumptions  of  Lauchle  (1981)  and  Somette  and  Lagier 
(1984)  have  on  estimation  of  the  sound. 

3.  The  effect  of  spot  intnaction,  which  occurs  in  an  actual  natural 
transitional  flow,  on  the  sound  produced  by  an  individual  turbuloit  ^t. 


16 


In  order  to  answer  these  questions  the  following  investigation  has  been 
performed.  A  turbulent  spot  was  gen^ted  in  a  laminar  flat  plate  boundary  layo*  by 
a  vigorous  impulsive  mass  injection  from  a  small  hole.  Measurements  were  taken  of 
the  velocity  field  fluctuation  caused  by  the  passage  of  the  spot  for  a  freestieam 
velocity  U  =  41  cm/s.  Specifically,  measurements  were  made  of  the  unsteady 
streamwise  velocity  profiles,  both  on  and  off  the  spot  coiterline,  at  five  streamwise 
locations  in  the  boundary  layer.  From  these  measurements  the  boundaries  of  the 
turbulent  spot  and  its  growth  rates  were  determined.  Using  conservation  of  mass  the 
local  rate  of  change  of  mass  flux  ova*  the  q>an  was  also  deduced.  The  rise  time,  tj, 
was  then  determined  using  the  local  rate  of  change  of  the  streamwise  mass  flux.  The 
normal  velocity  at  a  location  just  above  the  maximum  qx)t  height  was  also  measured 
directly  for  comparison  with  the  integral  analysis.  This  data  set,  then,  provides 
information  concerning  the  ensemble-averaged  unsteady  velocity  profile  over  a 
significant  portion  of  the  developmoit  of  the  turbulent  spot,  along  with  the  map  of  the 
potential  flow  perturbation  velocity  normal  to  the  wall,  for  one  set  of  conditions  in 
one  facility.  Such  a  complete  and  consistent  description  of  the  flow  quantities  needed 
for  the  Liepmann  acoustic  analogy  is  not  available  in  any  previous  experimental 
investigation  of  turbulent  spots. 

The  effect  of  turbuloit  spot  interaction  on  the  spot  growth  rates  was 
determined  next.  This  was  accomplished  by  goimting  a  second  spot  at  a  given  delay 
time  after  the  first.  The  streamwise  profiles  on  the  coiterline  were  measured  in  order 
to  deduce  the  same  streamwise  growth  parameters  as  in  the  baseline  case.  In  addition, 
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the  normal  velocity  was  measured  to  obtain  the  normal  mass  flux  distribution.  These 
measurements,  combined  with  the  single  spot  case,  allowed  an  assessment  of  the 
attenuation  of  the  growth  rates  of  a  spot  due  to  the  presence  of  another  spot. 

Finally,  the  sound  radiation  from  a  turbuloit  spot  was  estimated.  The  behavior 
of  the  normal  velocity  seen  in  the  measuremoits  enabled  the  normal  mass  flux 
distribution  of  a  turbulent  spot  to  be  described  by  a  superposition  of  simple  sources. 
This  arrangement  in  turn  led  to  a  scaling  analysis  of  the  turbuloit  spot  sound 
radiation.  This  analysis  was  used  in  conjunction  with  the  experimoital  data  to 
estimate  how  the  spot  generates  sound  as  it  grows.  Using  some  simple  assumptions. 


these  results  were  then  extmded  to  the  natural  transition  case. 


CHAFTER2 


ANALYTICAL  FOUNDATIONS 


In  order  to  estimate  the  sound  radiated  from  a  boundary  layar  flow,  it  is 
necessary  to  relate  the  farfield  pressure  radiation  to  motions  whidi  occur  in  the 
boundary  layer.  Information  about  the  unsteady  fluid  motion  in  the  boundary  layer 
may  be  dmived  either  from  calculation  mr  from  experimental  data.  The  rdationship 
b^ween  the  source  and  its  effect  on  the  farfield  has  been  presented  in  two  ways,  the 
{qrproach  of  Ughthill  (1952)  and  the  ry>proach  of  Liq>mann  (1954).  Lighthill  (1952) 
treated  the  seomd  derivatives  of  die  fluctuating  Reynolds  stresses  as  quadrupole 
sources  in  the  wave  equation  for  the  acoustic  radiation  (see  ntillips  (1955),  Powell 
(1960),  Hardin  (1991)).  This  information  is  not  easily  accessible  experimentally 
because  it  requires  measurement  of  the  sectmd  doivatives  of  second  and  higher  order 
fluctuating  vdodty  products  simultaneously  over  die  entire  source  field. 

An  alternative,  but  seldom  used,  aiqnoach  is  that  of  Liqimann  (1954),  which 
treats  the  potential  flow  fluctuations  direcdy  above  die  boundary  layer  as  a  boundary 
condititm  to  the  wave  equation.  A  similar  tq^iroach  was  later  used  by  Tam  (1975)  to 
analyze  the  sound  radiated  from  a  turbulent  boundary  layer.  Potential  flow 
fluctuations  may  be  dorived  from  fluctuations  in  the  boundary  layer  di^lacement 
thickness.  Using  concqits  from  slender  body  theory,  a  boundary  condition  for  die 
acoustic  wave  equation  may  be  constructed  from  the  known  unsteady  displacement 
thickness  bdiavior.  It  will  be  shown  in  this  chapter  that  while  the  form  for  the 
farfield  pressure  is  conveniendy  writtra  in  tmms  of  the  diiqilacement  thickness. 
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calculation  of  the  displacement  thickness  from  three-dimensional,  unsteady  velocity 
field  data  is  not  straightforward.  This  problem  has  motivated  the  use  of  a  related 
quantity,  the  mass  flux  deficit.  In  this  work,  the  unsteady  boundary  condition  for  the 
wave  equation  is  constructed  using  the  unsteady  mass  flux  deficit  obtained  from  the 
experimental  data. 

For  a  transitional  boundary  layo*  the  intermittent  nature  of  the  source  flow 
presents  a  situation  where  the  Liq)mann  ^roach  has  advantages  over  that  of 
lighthill  for  predicting  flow  noise.  The  distinct  changes  in  the  fiowfield  which  takt 
place  during  the  passage  of  a  turbulent  spot  allow  the  large-scale  fluctuations  of  the 
displacemmt  thickness  to  be  expressed  in  a  simple  form.  This  form  of  the  unsteady 
displacement  thickness  includes  parameters  which  may  be  derived  from  experimental 
data. 

2.1  THE  ACOUSTIC  BOUNDARY  VALUE  PROBLEM 

The  physical  system  investigated  here  was  a  uniformly  moving,  viscous  fluid 
flowing  over  a  flat  plate  with  a  uniform  velocity  U  ^  c,  where  c  is  the  ^)eed  of 
sound  for  the  fluid  which  is  uniform  throughout  (see  Figure  2.  la).  The  origin  of  the 
coordinate  system  is  at  the  leading  edge  of  the  flat  plate.  In  the  region  of  the  fluid 
near  to  the  wall  a  laminar  boundary  layer  is  formed.  A  turbuloit  spot  is  goimted 
by  a  point  impulsive  disturbance  at  a  location  x=Xo,  y=‘0,  z=0  (z=0  is  subsequently 
referred  to  as  the  centerline).  The  acoustic  source  n^on  is  the  portion  of  the 
boundary  layer  which  is  unsteady  due  to  the  presence  of  the  turbulent  spot.  For 
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purposes  of  analysis  the  fluid  region  is  divided  into  two  portions,  denoted  V  and  V«, 
as  shown  in  Figures  2.1a  and  2.1b.  The  inviscid  part  of  the  flow,  denoted  by  V,  is 
bounded  on  one  side  by  the  hemispherical  surface  S  located  at  a  radius  R^d,  and  on 
the  other  by  the  surface  S«,  which  is  located  outside  the  boundary  layer,  i.e., 
y>5(x,z,t),  where  8  is  the  boundary  layer  thickness.  The  viscous  flow  region,  V,, 
is  the  boundary  layer  itself,  bounded  by  the  surfaces  y  =  0  and  Sg.  The  radiation  of 
sound  from  the  boundary  layer  is  det^mined  by  relating  the  "nearfield”  disturbances 
on  S,  to  a  farfield  point  in  V  through  the  wave  equation.  Exp^moital  measuremrats 
inside  V,  provide  the  information  for  constructing  the  boundary  condition  on  S,.  The 
positions  of  an  observer  in  V  and  a  source  point  on  are  d«ioted  by  x  and  x', 
respectively,  while  the  distance  from  a  source  point  to  an  observer,  r  =  |  x'-x  | ,  is  the 
magnitude  of  the  vector  r  =  x'-x,  as  shown  in  Figure  2.1b. 

For  an  observer  anywhere  in  V,  but  far  enough  away  from  the  boundary 
layer,  the  sound  generated  by  the  tuibuloit  spot  is  described  by  the  acoustic  wave 
equation  which  has  been  linearized  about  the  steady  flow  : 


^  -  c^V^^  =0  in  V 
at^ 


(2.1) 
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(a)  Geometry  of  the  acoustic  boundary  value  problem,  showing  the  definiticms 
ofV,  V,. 


(b)  Detail  of  geometry  of  acoustic  boundary  valuee  problem  showing  the  r^cm  near 
the  source,  applied  to  a  transitional  bounds^  layer  flow. 


Figure  2.1  Geometry  of  the  acoustic  boundary  value  problem. 
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where  4>  is  the  acoustic  potratial.  The  boundary  condition  for  the  acoustic  potential 
is  given  by: 

^  -  v.(x,z,t)  on  S,  (2.2) 

where  v,  is  the  potential  flow  perturbation  velocity  due  to  the  presrace  of  the 
boundary  layer.  In  addition,  4>  satisfies  the  Sommerfeld  radiation  condition  on  S. 
Note  that  the  convected  wave  equation  is  not  bdng  used  since  M  <  1. 

Liq)niann  (1954)  used  this  approach  to  describe  the  sound  radiation  from  a 
turbulent  boundary  layer.  He  began  his  analysis  with  the  solution  to  the  boundary 
value  problem  just  described  : 

f  f  iZJ) 

2Tj,J  c  |x-x'| 

where  dS'  =  dx'dz'  is  the  differoitial  area  on  Sq.  For  the  inviscid,  moving  portion 
of  the  flow  in  V ,  the  perturbation  pressure  in  the  freestream  is  given  by  the  Bernoulli 
equation  linearized  about  the  steady  flow  condition: 

p(ii.t)  -  *  u4i]  a.4) 

dt  dx 


Substituting  Equation  2.3  into  this  relation,  he  obtained  the  expression  for  a  subsonic 
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freestream  over  the  planar  suifu%  S«  on  which  there  is  a  prescribed  normal  velocity 
distribution: 


HQ/  d  V 

_ _ {[ _ t]  -  U  [V  1 

|x-x'|  at  "  |x-x'|2 


} 


(2.5) 


where  the  notation  [  ]  indicates  a  quantity  evaluated  at  the  source  location  x'  and  at 
the  appropriate  retarded  time  t  -  |x'  -  xj/c.  At  lai;ge  r  =  |x'  -  x| ,  the  second  term 
is  n^ligible  compared  to  the  first,  so  that  it  may  be  dropped  whoi  considoing  the 
farfield  behavior.  This  reduces  the  expression  for  the  farfield  pressure  to: 


P(x,0 


dS^  .av 
|x-x'|  at 


(2.5) 


It  should  be  pointed  out  that  the  second  term  in  Equation  2.5  divided  by  -pU  is 
idoitical  to  that  used  by  Van  Atta  er  a/.(1982)  to  estimate  the  fluctuating  potential 
flow  velocity  normal  to  the  wall  at  a  hdght  just  above  the  boundary  layer. 

If  we  move  the  density  inside  the  int^ral  in  Equation  2.6,  we  may  interpret 
the  resulting  equation  as  the  time  derivative  of  the  unsteady  mass  flux  through  for 
each  point  on  S«  which  contributes  to  v.,  evaluated  at  the  appropriate  retarded  time. 
The  radiaticm  of  sound  from  the  boundary  layo-  is  then  d^ermined  by  the  form  which 
V.  takes.  If  v,  integrated  over  S  changes  with  time,  thra  the  flow  acts  as  a  monqwle 
type  source.  If,  however,  the  int^rated  mass  flux  is  i^tically  xero,  then  the  flow 
acts  as  a  higher  order  source  such  as  a  quadrupole.  The  quantitative  behavior  of  v. 
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may  be  found  experimentally  either  by  direct  measurement  or  by  measuremoits  of  the 
unsteady  velocity  in  the  boundary  layer  itself,  making  use  of  the  continuity  equation. 
To  accomplish  the  latter,  it  is  necessary  to  consider  more  closely  the  relationship 
between  the  unsteady  displacement  thickness  and  the  potoitial  flow  disturbance. 

2.2  BOUNDARY  CONDITION  FOR  ACOUSTIC  EQUATION 

By  measuring  the  flow  inside  the  boundary  layer,  it  is  possible  to  determine 
the  bdiavior  of  the  flow  on  So-  As  stated  in  Equation  2.2,  the  distribution  of  velocity 
normal  to  S,  constitutes  a  boundary  condition  for  the  acoustic  potential.  This  is 
accomplished  in  one  of  two  ways:  First,  the  velocity  component  normal  to  the  wall 
outside  the  boundary  layer  may  be  measured  directly.  Second,  the  instantaneous 
displacement  thickness,  d*(x,z,t),  may  be  calculated  at  every  point  for  every  instant 
using  the  experimental  data.  By  modeling  the  effect  of  the  boundary  layer  on  the 
outer  flow  as  that  of  an  additional  wall  thickness  y  =  5*(x,z,t),  it  is  possible  to 
position  S«  on  this  surface.  Thoi,  invoking  concepts  from  sloider  body  theory,  the 
boundary  condition  may  be  "transferred”  to  the  wall  (y=0). 

2.2.1  The  Dfeplacement  Thickness 

The  displacement  thickness.  S',  is  a  quantity  that  expresses  the  amount  of 
streamwise  mass  flux  that  is  lost  to  a  uniform  freestream  due  to  the  viscous  action  in 
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the  boundary  layer.  The  simplest  exampie  is  a  two-dimensional,  steady  laminar  flat 
plate  boundary  layer.  In  this  case  the  displacement  thickness  is  given  by: 

fix)  .  f(l  -  Ji^xiy  (i-T) 

At  any  stieamwise  location  in  a  flat  plate  boundary  layer,  the  displacement  thickness 
is  the  height  that  a  streamline  is  displaced  from  its  hdght  H  >  6  upstream  of  the  plate. 
From  the  point  of  view  of  the  potratial  flow,  the  displacement  effect  of  the  boundary 
layer  is  equivalent  to  a  uniform  inviscid  flow  slij^g  over  a  body  whose  shape  is 
adjusted  for  the  displacement  thickness.  For  a  three-dimensional  boundary  layer,  the 
di^lacement  thickness  is  a  function  of  changes  in  both  the  streamwise  and  the 
spanwise  mass  flux  in  the  boundary  layer.  These  statemcmts  follow  directly  from  the 
principle  of  conservation  of  mass  for  an  incompressible  flow.  For  example,  a 
decrease  in  the  stieamwise  mass  flux  wiU  be  reflected  in  an  increase  in  the  mass  flux 
normal  to  the  plate  (for  a  two-dimoisitmal  flow).  The  product  of  the  displacement 
thickness  and  the  fineestream  velocity  exactly  equals  the  streamwise  mass  flux  deficit, 
which  is  given  by: 

0» 

mass  flux  deficit  -  p  |  (U  -  u(x,y))dy  (2>8) 

for  a  steady,  two-dimensional  boundary  layer.  If  we  extend  this  description  of 
continuity  to  the  unsteady  flow  in  dther  tranutional  or  turbulent  boundary  layers,  we 
can  see  that  fluctuations  in  the  unsteady  displacement  thickness  reflect  the  exchange 
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of  mass  flux  between  fluctuations  in  flow  parallel  to  the  body  and  fluctuations  in  flow 
normal  to  the  body.  These  fluctuations  in  normal  velocity  are  to  be  used  as  the 
boundary  condition  for  the  acoustic  potential. 

2.2.2  Cakulating  Normal  Velocity  from  the  Displacement  Thickness 

While  the  displacemoit  thickness  is  a  useful  quantity  for  constructing  a 
boundary  condition  for  the  acoustic  equation,  it  is  not  necessarily  derived  from 
velocity  data  in  a  straightforward  way  as  Liq)nuum  indicated.  The  analytical  work  that 
has  proceeded  on  the  basis  established  by  liq>mann  that  is  relevant  to  this  work  has 
continued  to  use  the  displacement  thickness.  Because  of  the  difficulty  in  estimating 
the  displacement  thickness  using  data,  the  present  analysis  has  beat  based  on  a  related 
quantity,  the  mass  flux  deficit  (MFD).  This  quantity  relates  the  normal  vdocity 
perturbation  to  the  unsteady  boundary  layer  data  in  a  direct  manner,  while  the 
displacemrat  thickness  does  not,  as  will  be  shown  in  this  section. 

Expressing  the  fluctuating  normal  velocity  in  terms  of  an  unsteady 
displacement  thickness  proceeds  according  to  the  perturbation  theory  for  slender 
bodies  (Liepmann  &  Roshko,  1957).  Within  the  framework  of  this  theory,  the 
component  of  velocity  normal  to  the  freestream  flow  direction  of  the  flow  over  a  body 
is  the  result  of  the  body  shtqre  perturbing  the  invisdd,  uniform  flow.  Liqrmaim  & 
Roshko  (1957)  modeled  this  displacement  effect  as  an  effective  piston  motion  imposed 
on  steady,  potential  flows.  In  their  discussion  of  the  flow  of  visctHis  fluids,  they  note 
that  the  theory  can  be  extended  by  simply  changing  the  body  shsqre  by  the  addition  of 
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the  displacement  thickness,  assuming  that  there  is  no  boundary  layer  separation  on  the 
body.  For  a  steady,  two-dimensional  boundary  layer,  the  normal  perturbation  velocity 
is  entirely  due  to  the  effect  of  the  displacement  thickness  and  is  given  by: 

*(X) 

V.(x)|,^  -  .  uA  I  (1  -  2^)  dy  0-») 


where  x  is  the  streamwise  coordinate,  U  is  the  freestream  velocity,  d(x)  is  the 
boundary  layer  thickness,  and  5*(x)  is  the  displacement  thickness.  Note  that  in  the 
usual  definition  of  displacement  thickness  the  upper  limit  on  the  integral  is  infinity, 
while  here  the  integral  only  goes  up  to  the  boundary  layer  height.  Since  the 
displacement  thickness  is  defined  in  tarms  of  the  amount  of  mass  flux  lost  to  the 
incoming  stream  due  to  motion  in  the  boundary  layer,  the  upper  limit  on  the  integral 
is,  strictly  speaking,  the  boundary  layer  height.  In  a  steady,  two-dimoisional  flow, 
for  which  this  definition  is  illustrated,  the  nature  of  the  flow  allows  the  limit  of 
integration  to  be  extended  to  infinity  because  the  mass  flux  deficit  in  the  freestream 
is  entirely  negligible.  While  the  definition  of  5*  given  in  Equation  2.9  holds  for 
situations  other  than  a  two-dimensional,  steady  boundary  layer  flow,  the  relationship 
between  d*  and  the  mass  flux  deficit  is  not  always  so  straightforward. 

For  an  unsteady  two-dimoisional  boundary  layer  with  a  uniform,  steady 
potential  flow  upstream  of  the  plate,  the  equation  describing  the  normal  velocity  also 
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includes  an  unsteady  teim: 


«D 


dx 


(2.10) 


The  first  equation  was  given  by  Laufer  er  a/.(1964)  and  results  from  kineniatic  theory 
by  enforcing  the  condition  that  the  normal  velocity  of  the  fluid  relative  to  the  effective 
"wall"  (y  =  d*(x,t))  is  zero.  The  vdocity  of  the  fluid  is  tangent  to  the  moving  "wall" 
at  any  instant: 

^  ♦  (U  ‘Vy* «  0  (2.11) 

dt 

wh«e  f  =  d*(x,t)  >  y  is  the  surface  describing  the  effective  body  sluqw  due  to  the 
displacemoit  effect  and  U  is  the  vdodty  vector  U  =  (  U+u,  v,  w)'’'.  The  potential 
flow  poturbation  velodties  u,  v,  and  w  due  to  the  flow  over  the  body  are  snudl  such 
that: 


llUL  <  I 

U’U’U 


(2.12) 


The  resulting  unsteady  term  exhibits  the  motion  of  die  flow  streamlines.  The  second 


of  Equations  2.10  has  been  derived  by  intending  the  continuity  equation  for  an 
incompressible  flow.  The  relationship  bdween  the  two  iqiproadies  is  complicated  in 
two  ways.  First,  the  normal  vdodty  from  the  viscous/invisdd  intnaction  approach 
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is  given  for  y=5*,  while  the  normal  velocity  obtained  from  int^rating  the  continuity 
equation  is  evaluated  at  the  edge  of  the  boundary  layer.  Second,  the  ready 
identification  of  d*U  with  the  streamwise  derivative  of  the  mass  flux  deficit  has  bera 
lost.  However,  at  any  instant  the  normal  velocity  at  the  edge  of  the  boundary  layer 
is  still  given  by  the  streamwise  giadirat  of  the  streamwise  mass  flux  deficit. 

The  situation  is  complicated  further  in  a  three  dimensional  boundary  layer 
where  the  fireestream  is  uniform  and  steady  upstream  of  the  plate.  In  this  case  the 
normal  velocity  is  also  affected  by  spanwise  flow: 


v.(x.J,t)|,^  .  U{  ±1  (1  -  * 

<00 

_  a  f  W(x,y,z,t)^„  I 


(2.13) 


This  relation  was  givra  by  Lighthill  (1958).  The  three-dimensional  analog  to  the  first 
of  Equations  2. 10  is  given  by: 

v.(x.S-.z.t)  .  0.14) 

“  at  ax  dz 

This  expression  includes  the  contribution  of  the  spanwise  mass  flux  to  the  normal 
velocity. 


From  this  discussion  it  is  a^)parait  that  deducing  the  unsteady  normal  velocity 
perturbation  from  unsteady  boundary  layer  velocity  measurements  is  unnecessarily 
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complicated  by  referaice  to  the  displacement  thickness.  Instead,  it  is  more 
straightforward  to  proceed  from  velocity  profile  measurements  to  the  normal  velocity 
perturbation  through  integrating  the  continuity  equation.  For  this  reason  the  mass  flux 
deficit  is  used  here  instead  of  the  displacemoit  thickness,  as  in  Equation  2.13.  The 
use  of  the  displacemoit  thickness  may  work  well  for  a  formal  derivation  but  is  not 
well  suited  for  direct  sqtplication  to  experimental  data. 


2.2.3  Transfer  of  the  Boundary  Condition  to  the  WaU 

The  boundary  condition  for  the  acoustic  equation  (Equation  2.2)  may  be 
rewritten  in  terms  of  the  velocity  field  in  the  boundary  layer: 


^(x,8,z,t)  »  v,(x,«,z,t)  » 

ay 

.U^j(l-if^)dy. 

«  * 

-  -^1  w(x,y,z,t)dy 


(2.15) 


for  $•  defined  as  the  surfiu»  y=a(x,z,t).  For  simplicity  in  the  analysis,  the  boundary 
condition  may  be  transferred  to  y =0  using  ideas  from  slender-body  theory  (Liepmann 
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and  Roshko,  1957).  This  is  accomplished  by  expanding  Vo(x,6,z,t)  in  a  Taylor  series 
in  y  about  y=0: 

dv 

v,(x,8,z,t)  =  v.(x,0,2,t)  +  (•^)ly.o5 
By  neglecting  terms  after  the  first,  it  is  possible  to  state: 

4^(x, 0,1.1)  ”  v.(x,0,z,t)  - 

.  0.17) 

In  slendo’  body  theory,  this  stq)  is  takmi  before  the  analysis  is  carried  out  but  is 
confirmed  after  the  solution  is  in  hand.  This  involves  showing  that  limy.^dv/dy>-^), 
as  shown  by  Karamcheti  (1966).  Thus  far  a  formal  solution  for  the  boundary  layer 
noise  problem  has  been  developed  using  the  Liq)mann  analogy.  The  next  stq)  is  to 
supply  the  detailed  behavior  of  the  boundary  layer  velocity  field  so  that  the  boundary 
condition  2. 17  may  be  constructed  and  the  acoustic  equation  solved. 


2.3  BOUNDARY  CONDITION  FOR  A  TURBULENT  SPOT 

In  this  section,  the  formal  expressions  for  the  radiated  sound  which  motivate 
the  boundary  layer  measuremoits  are  presented.  The  details  of  the  sound  estimation 
based  on  these  expressions  using  the  data  are  presented  in  Cluq>ter  5. 

Lauchle  (1981)  analyzed  the  sound  production  from  the  transition  r^on  of  a 
flat  plate  boundary  layer  by  modeling  the  displacement  thickness  fluctuatitms.  Using 
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the  sq>proach  of  Libby  (1975),  Lauchle  (1981)  wrote  a  description  of  the  instantaneous 
displacement  thickness  as  a  weighted  sum: 

5‘(x,t)  =  I(x,t)6’(x)  +  (1  -  I(x,t))6  jx)  <2.18) 

where  I(x,t)  is  the  indicator  fiincticm.  Hiis  quantity  is  defined  as  follows:  it  takes 
two  values,  zero  when  the  boundary  layer  at  the  location  x  is  laminar  at  time  t,  and 
one  when  it  is  turbulent.  Also,  the  time  avnage  of  the  indicator  fimctimi  is  the  equal 
to  the  intermittency  : 

ITxJJ  =  7(x)  <2*1^ 

From  this  definition  it  can  be  seoi  that  the  time  average  of  Equation  2. 17  yields  the 
intermittotcy-weighted  form  of  the  displacement  thickness  established  by  Dhawan  and 
Naiasimha  (1958): 


l\x)  -  y(n)S^  *  (1  -  7(x))j' 

This  indicator  function-weighting  approach  assumes  that  the  displacemoit  thickness 
has  the  value  it  would  have  in  a  turbulent  boundary  layer  in  the  interior  of  a  turbulent 
spot,  and  the  value  it  would  have  in  the  undisturbed  laminar  boundary  layer  outside 
the  spot.  For  any  location  x  in  the  boundary  layer,  then,  when  a  turbulent  spot  is 
over  that  location,  d.(x,t)  >=  (d*)T,  othowise,  5*(x,t)  »  (5\.  The  instantaneous 
displacemoit  thickness  may  be  rewritten  as: 


where  A5*(x)  is  the  difference  b^een  the  values  of  the  laminar  and  turbuloit 
displacement  thicknesses  at  a  givoi  location  x.  Lauchle  further  assumed  that  the 
transition  between  the  laminar  and  turbuloit  states  was  not  instantaneous,  but  consisted 
of  an  accommodation  region  at  the  edges  of  the  spot  which  followed  an  exponratial 
evolution,  characterized  by  the  rise  tinm  as  shown  in  Figure  2.2.  Spatial  variations 
in  tj  were  not  considered.  The  unsteady  normal  velocity  was  related  to  fluctuations 
in  5*  through: 


da’ 

V.(x.t)  -  ^  - 
.  A«'(x)£(M) 


0.22) 


after  Lauf^,  et  al.,  (1964),  who  considered  the  sound  radiated  from  a  turbulent 
boundary  layer  to  an  obsover  movmg  at  U.  Substituting  this  expressicm  into 
Equation  2.6,  we  obtain  for  the  sound  radiation; 


P<M)  -  f  [  A»‘(x')t.^l  003) 

2t  ^  J  |x'  -  x|  dt* 

This  result  was  also  obtained  by  Audet  ct  a/.(1989)  from  a  Ughthill  analogy  aiq[>roach 
where  they  applied  indicator  functitm  wdghting  to  the  Ughthill  stress  tensor. 
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using  assumptions  similar  to  those  of  Lauchle  (1981).  It  should  be  mentioned  that 
Lauchle’s  (1981)  result  was  different  than  this  result  in  that  both  the  time  derivative 
and  the  (1/r)  term  were  outside  the  integral.  These  differraces  are  consistent  with  the 
assumptions  of  source  compactness  made  by  Lauchle  (1981). 

Thus,  using  an  indicator  function-weighted  description  for  the  displacement 
thickness  allows  the  sound  radiation  problem  to  be  formulated  in  terms  of  a  rise  time, 
(tj),  a  function  describing  the  change  in  displacement  thickness  at  each  streamwise 
location,  A5*(x)  (Lauchle  (1981)  neglected  the  variation  in  A5*  because  of  the  small 
streamwise  extent  of  the  transition  zone),  and  the  space-time  behavior  of  the  indicator 
function.  Since  Lauchle  (1981)  was  interested  in  describing  the  acoustic  emissions  of 
the  oitire  transition  region,  he  constructed  a  space-time  correlation  function  for  I(x,t) 
to  describe  its  behavior  over  the  extent  of  the  transition  zone.  A  description  of  the 
indicator  function  for  an  isolated  spot  may  be  expressed  in  a  more  detmministic 
manner  using  measurements  of  the  arrival  times  of  the  laminar/turbulrat  interfaces. 
Therefore,  to  implement  this  model  using  experimental  data,  it  is  necessary  to 
measure  Ad*(x),  tj,  and  the  space-time  behavior  of  I(x,t).  This  has  been  accomplished 
through  measurements  of  the  mass  flux  deficit  from  which  these  quantities  have  been 
calculated  (see  Ctuq)ter  4).  The  estimation  of  the  sound  radiation  from  transitional 
flows  through  the  Liepmann  analogy  using  the  experimental  data  is  presented  in 
Chapter  5. 


CHAPTERS 


EXPERIMENTAL  APPARATUS 


This  section  details  the  2q)paratus  used  to  carry  out  the  experimoits.  The 
design  and  construction  of  the  Laminar  Flow  Water  Channel  and  the  measurement 
system  are  described.  The  performance  of  the  facility  and  its  components  is  discussed 
along  with  some  suggested  improvements. 

3.1  WATER  CHANNEL  FACILITY 

All  of  the  experiments  for  this  study  took  place  in  the  Am^nce  Engineering 
Laminar  Flow  Water  Channel,  a  new  facility  constructed  over  a  period  of  two  and  a 
half  years  for  these  experimoits.  The  facility  is  located  in  the  DqKutmait  of 
Aerospace  Engineering  at  Penn  State.  Tl»  flow  loop  is  a  closed  circuit  which  holds 
sqqnoximatdy  3500  gallons.  All  but  cme  of  the  experiments  were  carried  out  at  a  test 
section  flow  speed  of  0.41  m/s,  although  the  maximum  speed  is  much  higho’.  The 
important  features  of  the  flow  loop  may  be  seen  in  Figure  3.1.  In  this  diagram  the 
flow  moves  in  a  clockwise  direction  around  the  loop.  The  first  important  feature  is 
the  test  section  where  the  experiments  took  place.  Downstream  of  the  test  section  is 
the  pump.  Continuing  in  the  flow  direction,  there  is  a  return  1^  consisting  of  a 
system  of  diffusa^,  a  turnaround  duct,  and  a  wide-angle  diffesa-.  The  return  1^ 
decelerates  the  flow  so  that  the  turbulence  management  causes  miiumal  flow  losses. 
Following  the  turbulence  management  thme  is  the  contraction  leading  into  the  test 


sectimi. 


2.00m 
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Figure  3.1  Wato*  channel  faK;ility.  Flow  moves  clockwise  around  the  loop. 
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3.1.1  Test  Section 

3.1.1.1  Floor  and  Walls 

The  test  section  is  3.66  m  long  and  has  a  nominal  cross  sectional  area  of  0.387 
m^,  which  varies  with  water  hdght  in  the  test  section.  The  test  section  bottom  is 
made  of  2.S4  cm  thick  hard  anodized  aluminum  and  the  walls  are  made  of  3  sectitms 
of  50.8  cm  X  1.22  m  X  1.27  cm  thick  glass  panes.  The  test  section  was  constructed 
on  a  welded  steel  frame  which  is  bolted  &st  to  the  laboratory  floor.  A  one-piece 
gasket  seals  the  joints  not  only  betweoi  the  base  plate  sections  but  also  the  glass  walls 
of  the  test  section.  For  access  to  the  test  section  from  boieath,  each  of  the  three  base 
plate  sections  has  a  27.9  cm  x  21.6  cm  plexiglass  window  insert  which  is  bolted  into 
the  base  plate  and  sealed  with  an  0-ring. 

On  either  side  of  the  test  section  is  welded  a  3.66  m  x  0. 102  m  C-beam,  which 
supports  a  2.54  cm  diameter  rail.  These  rails  support  an  instrummtation  rack  which 
is  constructed  of  10.2  cm  x  4.4  cm  sted  C-beam  sections.  The  stand  travds  alcmg 
the  rails  cm  sets  of  bearings.  Movemmit  altmg  the  rails  may  be  restricted  by 
tightening  a  bolt  in  the  bearing  mount  whidi  draws  the  bearing  mount  tight  around  the 
rail.  Loosening  these  bolts  allows  the  stand  to  move  fredy. 

The  test  section  is  connected  to  the  flow  loop  both  upstream  and  downstream 
by  aluminum  flanges.  These  flanges  are  attached  to  the  test  section  in  two  ways. 
Along  the  bottom  of  the  test  section,  the  flanges  are  bolted  to  the  base  plate,  while  tlm 
glass  side  walls  are  attached  (and  sealed)  to  the  flanges  using  siliccme  rubber  caulking. 
Rubber  gaskets  seal  the  joints  betweoi  the  flanges  and  their  adjacent  sections. 
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3.1.1.2  Flat  Plate  Test  Model 

The  test  plate  which  runs  the  length  of  the  test  section  is  constructed  from 
three  1.22  m  x  73.3  cm  sections  of  1.91  cm  aluminum  tool  and  jig  plate.  This  plate 
is  mounted  in  the  test  section  on  a  series  of  rails  at  a  height  of  12.7  cm  from  the  floor 
of  the  test  section.  This  positions  the  test  surface  14.6  cm  above  the  test  section  floor, 
leaving  as  much  as  31.1  cm  of  dq)th  beneath  the  free  surface.  Hie  support  rails 
beneath  the  surface  are  streamlined  on  their  leading  and  trailing  edges  to  help  reduce 
blockage  boieath  the  test  plate. 

A  12: 1  ellipse  is  machined  into  the  leading  edge  of  the  model.  In  addition,  a 
movable  tail  with  a  15.2  cm  chord  loigth  is  attached  to  the  trailing  edge  of  the  model 
to  enable  control  of  the  leading  edge  stagnation  point  location.  As  shown  in  Figure 
3.2,  two  dye  slots  are  mounted  in  the  model,  one  each  in  the  first  sections.  These 
slots  are  fed  from  a  manifold  which  is  connected  to  four  0.635  cm  copper  tubing  lines 
which  lead  to  a  head  tank  from  which  the  dye  is  dispensed.  The  dye  that  was  used 
in  these  experiments  consisted  of  a  Fluorescdn  Disodium  Salt  solution  which  emits 
green  light  when  excited  by  blue  light.  Also  machined  in  the  test  surface  are  5. 1  cm 
diameter  circular  ports  for  surface-mounted  sensors  and  three  20.3  cm  x  27.9  cm 
rectangular  windows  for  optical  access  from  below.  Plexiglass  inserts  are  machined 
and  mounted  into  these  ports  flush  to  the  test  surface.  The  rectangular  windows  in 
the  two  upstream  plate  sections  also  have  two  5. 1  cm  diameter  circular  ports  machined 
into  their  craterlines  and  plexiglass  inserts  are  available  for  these  ports  as  wdl. 
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Figure  3.2  Flat  plate  test  model. 
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A  provision  is  available  for  control  of  the  sidewall  contamination  of  the 
laminar  boundary  layer  by  adjustable  slots  at  the  side  edges  of  the  plate.  As  shown 
by  Motohashi  and  Blackwelder  (1983),  controlling  the  sidewall  contamination  should 
allow  for  a  greater  plate  area  exhibiting  laminar  flow.  This  is  accomplished  by 
causing  a  higher  pressure  on  the  upp^  side  of  the  plate  than  on  the  lower  side,  which 
allows  the  turbulent  boundary  layers  on  the  sidewalls  to  be  bled  away,  delaying 
destabilization  of  the  test  sur£ace  laminar  boundary  layer.  The  g^  between  the  edge 
of  the  plate  and  the  sidewall  is  nominally  1.27  cm  and  is  adjustable  down  to  zero  by 
using  four  nylon  tabs  on  each  side. 

3.1.1.3  First  and  Second  Comers 

The  test  section  is  followed  by  a  1.22  m  length  of  empty  channel  which  acts 
as  a  buffer  between  the  test  section  and  the  first  comer.  This  section  insures  that  any 
secondary  flow  induced  by  the  turn  will  have  a  reduced  influence  on  the  flow  in  the 
test  section.  Next  comes  the  first  comer.  The  flow  through  the  first  comer  is  aided 
in  turning  by  a  0.159  cm  thick  stainless  steel  perforated  plate  placed  at  a  45°  angle 
to  the  oncoming  flow.  This  plate  has  0.32  cm  poforations  and  a  64%  opoi  area 
ratio.  Beyond  this  plate  the  cross  sectional  area  is  then  reduced  through  the  second 
comer  to  that  of  the  pump  (856  cm^. 

3.1.2  Pump 

The  axial  flow  pump  utilizes  a  propeller  driven  by  a  7.5-HP  DC  electric 
motor.  A  schematic  of  the  pump  is  ^own  in  Figure  3.3.  The  power  of  the  motor 
is  transmitted  to  the  propeller  through  a  157  cm  long,  3.5  cm-diameter  stainless  steel 
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F^re  3.3  Detail  of  water  channel  pump. 
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shaft.  The  shaft  is  connected  to  the  motor  using  a  flexible  coupler.  The  propeller  is 
a  41  cm  brass  boat  propeller  which  has  bera  "tipped"  on  a  lathe  to  a  nominal  30.S  cm 
diameter  to  fit  in  the  30.S  cm  passage.  The  shaft  is  supported  in  two  places:  by  a  ball 
bearing  mounted  in  a  vibration-isolated  bracket  bolted  directly  to  the  floor,  and  by  a 
self-lubricating  graphite  bushing  located  in  the  flow  passage.  The  radial  loads  on 
these  supports  are  minimal,  and  the  outside  ball  bearing  supports  the  approximately 
45  N  thrust  load  from  the  propeller.  The  shaft  feedthrough,  where  the  shaft  passes 
through  the  wall  of  the  flow  passage,  consists  of  a  brass  tube  soldmed  to  a  brass  plate 
which  is  mounted  to  the  wall,  as  shown  in  Figure  3.4.  Sealing  this  wall  opening  was 
accomplished  by  covering  the  brass  plate  and  tube  with  fiberglass-reinforced  qx>xy 
resin.  On  the  outside  end  of  the  brass  tube  a  Chicago  Rawhide  CR-11171  seal  is 
mounted.  This  seal  has  worked  quite  effectively  in  sealing  this  potratially 
troublesome  joint.  To  prevent  motor  vibrations  from  transmitting  through  the 
laboratory  floor  into  the  water  channel,  the  pump  motor  is  mounted  on  a  heavy  steel 
base  which  in  turn  is  mounted  to  the  floor  using  VUct  CMB-7  vibration  isolators. 

3.1.3  Return  Leg 

The  return  leg  of  the  flow  loop  begins  immediately  downstream  of  the  pump. 
This  section  was  designed  to  increase  the  cross-sectional  area  of  the  flow  passage  from 
that  at  the  pump  to  that  of  the  settling  chamb^  as  efficiently  as  possible  while  Imqnng 
the  flow  in  the  test  section  steady.  The  design  was  constrained  by  the  geometry  of 
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Figure  3.4  Detail  showing  propdler  shaft  feedthrough  and  how  it  is  sealed. 
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both  the  room  where  the  facility  is  located  and  the  salvaged  sections  of  the  original 
settling  chamber. 

The  first  section  is  a  4.57  m  long  x  35.6  cm  high,  small-angle  (8**)  diffuser. 
This  diffuser  begins  at  a  joint  with  the  propeller  housing  which  is  made  of  two  flanged 
PVC  30.5  cm  pipes.  One  flange  is  connected  to  the  pump  and  one  to  the  diffusa:. 
The  joint  between  the  PVC  pipe  lengths  consists  of  a  pipe  clamp.  At  the  inlet  the 
diffuser  has  a  34.3  cm-diameter  circular  cross-section.  This  sh^  gradually  changes 
to  a  rectangular  cross  section  220  cm  downstream  of  the  joint.  The  outside 
construction  of  the  first  diffuser  is  a  rectangular  cross  section,  with  the  internal  shape 
formed  by  filling  the  comers  of  the  rectangular  duct  with  fiberglass-reinforced  qx>xy 
resin.  The  resin  was  used  to  seal  the  inside  surfiu:e  from  leaks  and  to  reinforce  the 
comer  joints  in  the  wooden  walls.  The  walls  are  also  reinforced  outside  the  flow 
passage  by  ribs  made  of  5.1  cm  x  1.7  cm  boards,  spaced  tq)proximately  every 
30.5  cm  in  the  axial  direction. 

The  diffuser  angle  and  length  woe  based  primarily  on  the  performance  data 
of  Japikse  (1984).  This  diffuser  is  not  within  the  transitory  stall  regime.  Transitory 
stall  would  be  best  for  pressure  recovery,  but  is  unaccq)table  from  a  steady  flow 
pospective.  An  analysis  was  carried  out  to  estimate  the  locations  of  any  flow 
sQMuations  in  the  diffuser  using  the  STAN6  boundary  layo*  code,  devdoped  by  W.M. 
Kays.  This  is  a  finite  difference  code  which  uses  a  mixing  length  turbulence  model. 
The  procedure  for  locating  flow  separation  was  as  follows:  First,  the  pressure 
gradient  given  by  a  uniform  invisdd  flow  through  the  diverging  duct  was  used  to 
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calculate  the  boundary  layer  up  to  sqxaration  (the  point  where  the  wall  shear  stress 
went  to  zero).  Then,  a  correction  for  the  wall  shape  based  on  the  calculated 
displacemoit  thickness  distribution  of  the  boundary  layer  was  used  to  recalculate  the 
pressure  gradient  distribution.  This  continued  iteratively  until  the  displacemrat 
thickness  matched  the  adjustment  to  the  duct  incorporated  in  the  code.  This  procedure 
did  not  account  for  any  three-dimensi<Hial  effects  due  to  variations  in  duct  cross 
sectional  sitapt  and  did  not  account  for  inlet  swirl  induced  by  the  propeller.  The 
swirl,  if  significant,  would  aid  in  delaying  sq>aration,  so  it  is  believed  that  the 
calculation  is  conservative  and  that  no  large-scale,  unsteady  sqjaradons  occur  in  the 
first  stage  of  the  diffuser.  Experimental  observation  has  not  revealed  any  large-scale 
unsteadiness  which  would  indicate  transitory  stall. 

The  second  stage  of  the  return  1^  is  a  wider  angle  diffuser,  as  shown  in 
Figure  3.5.  In  its  original  configuration,  the  diffuse  angle  was  bO**.  This  angle 
would  result  in  a  fully-separated  "jet”  flow,  so  the  angle  was  reduced  to  40**  by 
installing  wooden  inserts  to  change  the  stuq)e  of  the  flow  passage.  In  addition,  six 
80.6  cm  long  splitter  vanes  are  positioned  to  divide  the  duct  into  seven  passages. 
These  vanes  were  manufactured  from  marine-grade  plywood  and  coated  with 
fiberglass-reinforced  epoxy  resin.  The  leading  edges  are  rounded  to  prevent  leading 
edge  sqparation  and  the  trailing  edges  are  t2q)»ed  to  reduce  the  thickness  of  the 
trailing  wakes. 

The  positions  of  these  vanes  was  detmnined  using  the  criteria  set  by  Feil 
(1964).  In  goieral,  splitter  vanes  are  placed  such  that  they  divide  the  diffuser  passage 
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evenly.  In  doing  so,  the  vanes  lie  on  rays  emanating  from  the  point  of  intersection 
of  lines  drawn  extrading  the  di^ser  walls  into  the  inlet  duct  (point  O  in  Figure  3.6). 
Feil  advocates  placing  the  vanes  on  rays  emanating  from  a  "source-point"  O'  which 
is  upstream  of  point  O.  This  placement  helps  to  prevent  separation  at  the  location 
most  likely  to  sqrarate,  the  inlet  corner,  where  the  flow  curvature  is  the  highest. 
This  result  reflects  several  factors.  Fiiat,  the  vanes  help  to  turn  the  flow  around  the 
comer.  Second,  the  angle  is  smaller  in  the  outside  passages,  which  are  more  lilroly 
to  sqrarate  because  the  boundary  layer  on  the  sidewall  began  growing  far  upstream 
of  the  comer  where  the  vanes  begin.  In  this  way,  a  "crowding"  of  the  potential  flow 
in  the  di^ser  due  to  the  growing  boundary  layers  (especially  the  sidewalls)  helps  to 
retard  separation.  Lastly,  dividing  the  duct  into  smaller  passages  allows  any 
separations  that  do  occur  to  be  contained.  The  comras  at  the  diffuser  joint  were  also 
rounded  to  further  guard  against  separation. 

A  provision  for  boundary  layer  control  was  also  added  at  the  joint  ahead  of  the 
wide-angle  diffuser.  Here  two  slots  are  placed  in  the  wall.  The  filter  intake  was  a 
simple  slot,  as  shown  in  Figure  3.7.  Suction  at  this  location  removes  the  low 
momentum  fluid  near  the  wall  so  that  the  boundary  layer  is  less  likely  to  sqrarate  at 
the  comer.  Just  downstream  of  the  suction  slot  a  blowing  slot  (the  filtration  outlet) 
directs  a  jet  along  the  wall,  again,  to  re-raergize  the  boundary  layer.  A  filtration 
system  has  been  installed  to  process  the  water  through  this  loop. 

In  addition  to  the  vanes  just  mrationed,  two  rows  of  vanes  were  placed  just 
downstream  of  the  first  row.  These  vanes  help  to  turn  the  flow  through  the 
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Figure  3.5  Detail  of  the  compound  diffoser  showing  the  location  of  the  filtraticHi 
intake/outlet  and  the  splitter  vane  configuration. 


Figure  3.6  Diffuser  scnuce-point  di^lacement  (Feil,  1964).  Ideal  flow  in  diffuser 
looks  like  flow  from  source  at  O.  Vanes  set  on  rays  emanating  from  source  at  O' 
help  retard  sqnration  at  diffuser  joint  comers. 
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Figure  3.7  Schematic  of  filtration  inlet/outlet  used  as  boundary  layer  control. 
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decreasing  angle  of  the  flow  passage.  It  should  be  emphasized  that  the  design  of 
These  diffusers  optimizes  the  flow  passage  for  steadiness,  rather  than  for  pressure 
recovery.  Downstream  of  the  diffuser  the  flow  passes  through  two  90°  comers. 
The  flow  passage  at  this  location  is  2  m  wide  and  3S.6  cm  high.  The  comers  each 
have  a  15.2  cm  turning  radius.  To  help  prevrat  sqnration  of  the  flow  in  the  comers, 
six  PVC  turning  vanes  with  a  1.27  cm  nominal  thickness  and  a  21.6  cm  chord  length 
are  installed.  The  gs^s  between  the  vanes  ate  7. 1 1  cm,  so  that  the  chord-to-gs^  ratio 
is  f^roximately  3.  The  vanes  themselves  are  made  of  quarter  sections  of  a  30.5cm 
diameter  PVC  pipe  cut  loigthwise.  The  leading  edges  of  the  vanes  w^  rounded  and 
the  trailing  edges  were  t2q)ered  to  improve  the  flow  quality. 

The  final  diffuser  begins  at  the  exit  of  the  second  comer.  This  diffuser 
increases  the  flow  area  from  0.356  m^  to  1.14  m^  in  1.14  m  of  streamwise  distance. 
This  rapid  expansion  decelerates  the  flow  so  that  the  flow  conditioning  devices  in  the 
settling  section  (a  honeycomb  and  screws)  will  not  cause  a  large  head  loss  to  the 
flow.  The  diffuser  is  equipped  with  seven  perforated  plates  (0.159  cm  thickness, 
0.3175  cm  perforations,  64%  open  area)  which  act  to  re-oiergize  the  wall  boundary 
layers  and  prevent  them  from  sqKuating.  These  perforated  plates  also  act  to  break 
up  the  wakes  of  the  turning  vanes  so  that  the  flow  is  uniform  upon  leaving  the 
diffuser.  The  number  of  perforated  plat^  was  found  by  a  calculaticHi  similar  to  that 
done  for  the  lower  diffuser.  The  STAN6  code  was  run  until  separation,  and  then  the 
effect  of  a  perforated  plate  was  modeled  by  reducing  the  momentum  thickness 
Reynolds  number  by  50%  and  the  code  was  mn  iteratively  again  until  the  flow  was 
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found  to  separate.  The  separation  distance  was  found  to  increase  with  distance 
travelled  along  the  diffuser  axis.  The  choice  of  seven  plates  was  a  conservative 
estimate  of  plates  necessary  to  prevent  sq>aration  because  the  freestream  turbulence 
introduced  by  the  perforated  plates  also  aids  in  preventing  separation. 

Following  the  diffuser  is  the  settling  section  where  the  flow  conditioning  is 
installed.  First,  there  is  a  15.2  cm  Icmg  honeycomb,  which  has  0.32  cm-diameter 
cells.  The  plastic  honeycomb  was  manufactured  by  Plascore,  Inc.  Second,  there  is 
the  provision  for  fine-mesh  screens,  which  will  be  installed  at  a  future  date.  The  flow 
conditioning  is  completed  by  passing  the  flow  through  a  6:1  contraction  before  it 
enters  the  test  section. 

3.2  DATA  ACQUISITION 
3.2.1  Laser  Doppler  Velociinrter 

The  laser-Doppler  Velocimeter  used  for  all  measurements  was  a  Thermal 
Systems,  Inc.  (TSI)  one-component  system,  based  upon  a  15  mW  Helium-Neon  laser. 
A  schematic  of  the  system  components  is  shown  in  Figure  3.8.  The  laser  is  mounted 
vertically,  pointing  downwards,  on  the  travelling  instrumentation  stand.  A  right-angle 
mirror,  which  moves  with  the  optics,  directs  the  laser  beam  into  the  optical  train 
which  is  mounted  in  either  forward  or  backscatter  mode.  The  transmitting  q)tics  were 
mounted  so  that  the  beam  splitter  could  be  rotated  to  90*^  in  order  to  measure  the 
vertical  component  of  velocity  directly.  The  optical  axis  is  positioned  so  that  whoi 
measuring  the  streamwise  componoit  of  velocity  the  beams  approach  the  test  surface 
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Figure  3.8  Schematic  of  LDA  setup.  Laser  and  (:q)tics  are  mounted  to  the 
instrumentation  stand.  Flow  is  into  the  page. 
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at  a  slight  angle  and  are  not  blocked  whm  the  probe  volume  was  less  than  a  beam 
diameter  from  the  wall.  The  optical  train  is  been  mounted  to  a  Velmex  A4012 
spanwise  motion  traverse  that  has  15.24  cm  of  available  travel.  This  traverse  is 
driven  manually.  The  spanwise  travme  and  the  right-angle  mirror  that  directs  the 
laser  beam  into  the  LDV  optics  were  both  attached  to  an  aluminum  plate  which  was 
translated  vertically  using  a  Velmex  B4024  motorized  traverse.  The  vertical  traverse 
lead  screw  and  stepper  motor  provide  precise  motion  ~  1574  steps/cm.  This  motion 
is  controlled  by  a  Coitroid  Motion  Controller.  Commands  are  sent  to  the  Centroid 
controller  through  the  serial  port  of  the  lab  computer,  an  IBM-PC  AT.  A  Turbo 
Pascal  program,  which  also  controls  the  data  acquisition,  is  used  to  send  the 
commands.  The  combination  of  the  streamwise  travel  of  the  instrumoitation  stand, 
the  manual  spanwise  traverse  and  the  motorized  vertical  traverse  allows  three- 
dimensional  positioning  capability  for  the  LDA  system. 

3.2.2  Data  Acquisition 

The  signal  from  the  LDV  photodetector  was  fed  into  to  a  TSI  Model  1980B 
counter  processor,  equipped  with  a  Model  1988  Analog  Out  Module.  The  counter 
was  set  to  process  Doppler  bursts  with  at  least  16  fringe  detections.  The  time  series 
velocity  data  were  acquired  using  the  frequency-out  port,  which  delivmed  a  signal 
between  0  and  10  volts  proportional  to  the  binary  voltage  mantissa  of  the  last  validated 
Doppler  burst  processed  by  the  counter.  Thus,  the  analog  out  module  provided  a 
random  sample  and  hold  representation  of  the  velocity  time  series.  Using  the  sample- 
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and-hold  signal  automatically  corrects  for  velocity  bias  errors  in  the  average  velocity. 
The  analog  output  from  the  frequency-out  port  was  sampled  at  250  Hz  on  a  PC-based 
MetraByte  DAS- 16  analog-to-digital  conversion  board  cq)erating  in  the  polled  mode 
on  a  0-10  V  gain  setting. 

A  Turbo  Pascal  program  was  used  to  control  the  measurement  of  the  velocity 
profiles.  First,  the  program  asks  for  the  streamwise  and  qsanwise  locations  of  the 
probe  volume.  Thai,  it  asks  the  user  to  position  the  probe  volume  in  the  freestream 
so  that  the  freestream  velocity  at  that  location  can  be  measured.  Next,  the  counter 
settings  that  will  be  used  for  the  profile  measurement  are  entered  into  a  data  file. 
Once  the  freestream  velocity  has  been  measured,  the  usa  is  asked  to  position  the 
probe  at  the  wall  and  the  program  waits  until  the  user  tells  the  program  that  the  probe 
is  so  positioned.  The  user  positions  the  probe  by  telling  the  program  how  many  stq>s 
and  in  what  direction  to  move  the  traverse.  This  enables  the  user  to  position  the 
probe  at  the  wall  as  follows:  First,  a  shaip-cornered  block  of  known  height  is  placed 
on  the  test  plate.  Second,  the  probe  volume  is  brought  down  so  that  it  just  barely 
touches  the  corner  of  the  block.  Looking  at  the  oscilloscope  during  this  procedure 
allows  the  user  to  find  the  vertical  position  of  the  probe  volume  for  which  the 
maximum  amount  of  light  is  scattered  from  the  probe  volume.  This  ensures  a 
repeatable  means  of  finding  the  location  of  the  wall.  Once  the  wall  has  been  located, 
the  program  follows  the  procedure  given  in  the  previous  paragr^h  for  measuring  the 
velocity  profile  using  a  predefined  list  of  y-locations  for  measurement. 
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To  measure  the  unsteady  proves,  the  data  acquisition  must  be  triggered  based 
upon  the  generation  of  the  turbulent  spots.  Details  of  the  spot  generation  procedure 
are  given  in  the  next  section.  The  profiles  consisted  of  time  series  of  velocities 
measured  at  22  heights  above  the  plate.  The  distance  between  heights  was  scaled  on 
the  thickness  of  the  laminar  boundary  layer,  5i,,  at  the  location  in  question.  The  first 
12  measurement  heights  were  evenly  spaced,  separated  by  bj\l.  The  last  10  were 
also  evenly  spaced,  but  separated  by  bJZ.  Thus,  the  extent  of  the  vertical  traverse 
was  just  over  four  laminar  boundary  layer  heights.  At  each  measurement  location, 
53  spot  samples  were  taken.  This  number  was  found  to  be  sufficient  to  achieve  a 
convergence  of  an  ensemble  average.  Figure  3.9  shows  the  convergoice  of  the 
ensemble  average.  In  this  figure,  the  ordinate  is  the  time-average  of  the  absolute 
value  of  the  difference  between  the  average  of  n  samples  and  the  average  of  n-1 
samples.  Figure  3. 10  shows  the  oisemble-averaged  time  series  for  the  mean  and  the 
variance  of  the  streamwise  compon«itof  velocity  at  y=0.8  mm,  x=1.53  m,  z=0  m, 
and  U=40.6  cm/s.  Note  that  the  variance  outside  the  portion  of  the  time  series 
corresponding  to  the  turbulent  spot  is  low,  indicating  that  the  measurement  noise  is 
a  small  part  of  the  variance  associated  with  the  spot. 

The  process  of  generating  spots,  acquiring  the  time  series,  and  then  moving 
the  traverse  to  the  next  station  was  controlled  by  the  laboratory  computer.  The 
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Figure  3.9  Convergoice  of  the  oisemble  avoage. 
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for  generating  and  acquiring  velocity  time  series  data  from  a  single  tuibulmt  spot 
realization  was  as  follows: 

1.  Pulse  spot  generator 

2.  Acquire  analog  out  voltage  for  specified  time  (S-10  s)  at  250  Hz 

3.  Reset  spot  goierator  (To  oisuie  that  the  same  pump  motion  was 
followed  for  each  sample,  the  pump  was  reset  afrer  every  g«ieration 
event) 

4.  Allow  flow  to  settle  for  10  s  so  that  the  wake  of  the  previous  spot  clears 
the  measurement  location  in  order  to  prevoit  spot  interaction. 

Each  cycle  took  from  15  to  20  s. 

A  similar  sequence  of  evrats  was  used  to  measure  Vp,  the  potoitial  flow 
velocity  fluctuation  normal  to  the  wall.  For  this  measurement  the  transmitting  optics 
had  to  be  rotated  90**  so  that  the  probe  volume  fringes  were  oriented  parallel  to  the 
freestream.  A  Bragg  cell  was  used  to  provide  a  frequency  shift  of  100  kHz  in  order 
to  cause  the  interference  fringe  pattern  to  move  vertically.  This  technique  allowed  the 
fluctuation  of  the  normal  component  of  velocity  to  be  measured  directly.  The 
fluctuations  in  velocity  were  small  (typically  1%  of  the  freestream  velocity)  and 
individual  realizations  of  the  time  series  were  noisy.  The  peak  voltage  fluctuation 
registered  by  the  counter  corresponding  to  the  peak  velocity  amplitude  was  typically 
±0.09  volts  from  the  DC  value,  or  a  fluctuation  of  ±32  counts  in  the  A/D 
conversion.  Since  the  lowest  fluctuation  that  the  LV  should  measure  is  jq>proximately 
0.1%  of  the  freestream  velocity,  the  uncertainty  in  these  measuremoits  is 
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approximately  10%.  Improved  accuracy  was  obtained  by  subtracting  the  reading 
obtained  without  a  spot  present  to  eliminate  any  bias.  The  normal  velocity 
measurements  were  made  above  the  maximum  turbulmt  ^t  height,  at  approximately 
y  =  3.86l. 

3.3  SPOT  GENERATOR 

Turbulent  spots  were  generated  by  an  impulsive  jet  ejected  from  a  small  hole 
in  the  test  plate.  A  peristaltic  pump  driven  by  a  stepper  motor  produces  the  jet.  The 
pump  was  connected  by  0.635  cm  diameter  copper  tubing  to  the  disturbance  point,  a 
0.5  mm  hole  located  on  the  plate  coiterline  53  cm  from  the  test  surface  leading  edge. 
The  strength  of  the  pulse  may  be  varied  by  changing  the  motion  that  the  stepper  motor 
executes  in  driving  the  pump  to  produce  the  jet.  The  pulse  stroigth  for  the 
experimrats  was  determined  by  flow  visualization,  using  a  near*wall  sheet  of 
Fluorescein  dye  which  issued  from  the  most  upstream  dye  slot.  The  number  of  motor 
steps  for  a  disturbance  pulse  was  varied  and  the  resulting  turbulent  spots  observed. 
The  pulse  strength  was  chosen  by  determining  the  lowest  numb^  of  motor  stq)s 
needed  to  produce  a  spot  consistratly.  This  was  found  to  be  a  motor  motion  of  6 
steps,  or  0.03  revolutions.  The  impulsive  jet  introduced  an  additional  1.36  cc  of 
water  into  the  flow.  An  ensemble-averaged  time  trace  of  the  disturbance  signature  is 
shown  in  Figure  3.11.  The  disturbance  is  the  large,  sharp  peak,  which  is  followed 
by  smaller  oscillations.  These  oscillations  are  typical  of  an  impulsive  jet  gmerated 
by  a  square-wave  pulse  travelling  through  a  long  tube. 


3.4  WATER  CHANNEL  QUALIFICATION 
3.4.1  Water  Quality 

To  enable  repeatable  measurements  the  water  channel  must  be  kept  as  clean 
as  possible.  For  IDA  measuremeits  it  is  desirable  to  minimize  the  number  of 
particles  in  the  water  which  are  larger  than  a  few  microns  in  diameter.  Particulate 
matte  can  also  foul  the  honeycomb  and  any  turbulence  cmitrol  screens.  A  filtration 
system  was  installed  to  remove  large  particles  fiom  the  flow.  The  wate  was  removed 
from  the  main  flow  loop  through  the  boundary  laye  suction  slot  and  re-injected  just 
downstream,  as  shown  in  Figure  3.7.  It  was  passed  through  a  S-micron  filter  to 
remove  the  undesired  material.  Despite  this  measure,  a  great  deal  of  particulate 
matter,  largely  in  the  form  of  aluminum  corrosion,  was  always  present  in  the  flow. 
This  problem  motivated  periodic  fiu^ning  of  the  attire  system. 

The  corrosion  of  anodized  aluminum  surfaces  was  severe  in  some  cases. 
Nowhere  was  this  of  more  cmncan  than  the  test  surface.  Small  pinholes  in  the 
anodized  surface  allowed  the  test  plate  to  corrode  in  small  pits.  The  test  surface  had 
to  be  regularly  treated  to  slow  fiirther  corrosion.  This  treatmont  consisted  of  filling 
the  pits  with  Devcon  aluminum  putty,  sanding  the  putty  flush  to  the  surface  and  that 
applying  two  coats  of  car  wax  to  the  aluminum  surfaces.  The  wax  coats  lasted  up  to 
two  weeks  before  anotha  coat  was  needed.  This  treatment  was  found  to  be  effective 
in  protecting  the  test  surface. 
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3.4.2  Baseline  Flow 

LDA  surveys  revealed  the  quality  of  the  baseline  flow  in  die  test  sectiixi. 
Figure  3.12  shows  the  centerline  laminar  boundary  layer  profile  in  similarity 
coordinates  compared  to  the  Blasius  profile.  In  addition,  the  free^ream  velocity  was 
found  to  have  very  good  streamwise  and  spanwise  uniformity.  An  upper  bound  on 
the  turbulence  intensity  of  the  freestream  flow  was  found  to  be  0.3%  by  measuring 
the  variance  of  the  analog  signal  from  the  counter.  In  addition,  there  were  r^ons 
of  wall  contamination  which  spread  at  s^roximately  10**  from  the  wall,  limiting  the 
area  over  which  a  laminar  boundary  layer  existed  on  the  test  surface,  as  semi  in 
Figure  3.13.  These  r^ons,  first  noticed  by  Charters  (1943),  consisted  of  intermittent 
flow.  The  intmmittency  was  more  pronounced  at  lower  freestream  velocities.  The 
apex  of  the  triangle  formed  by  the  contamination  wedges  was  located  in  the  vicinity 
of  the  centerline.  This  point,  denoted  is  the  furthest  downstream  location  for 
which  uncontaminated  centerline  measurements  may  be  taken.  The  variation  in 
with  freestream  velocity  U  is  shown  in  Figure  3.14.  U  was  found  using  the  LDA 
and  x^  was  found  by  observing  the  most  upstream  location  of  the  disturbance  to  a 
sheet  of  dye  adjacent  to  the  wall  due  to  turbulmit  mixing. 

3.4.3  Suggested  Improvements 

The  recommended  improvemmits  to  the  water  channel  facility  may  be  divided 
into  two  categories.  These  are  the  quality  of  the  water  and  the  quality  of  the  flow. 
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Figure  3.13  Schematic  of  test  plate  flow  r^mes,  showing  the  sidewall 
contamination. 
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At  present,  experiments  in  the  water  channel  are  plagued  by  particulate  matter 
in  the  flow.  Although  a  filtration  system  has  been  installed,  it  has  not  been  entirely 
effective.  There  are  two  possible  reasons  for  this  failure:  1)  the  Alter  pump  is  not 
powerful  enough  to  produce  the  needed  suction  and  2)  the  access  to  the  particulate 
matter  is  limited  because  the  Altration  in/outlets  are  located  in  the  sidewalls  rather 
than  the  floor  of  the  lower  leg  difliiser.  To  remedy  the  first  possibility,  a  larger  pump 
should  be  installed.  Access  to  a  greater  proportion  of  the  flow  would  be  a  more 
ambitious  undertaking,  since  it  involves  modiflcations  to  the  flow  passage  structure. 
However,  a  more  permanent  (and  expulsive)  r^roach  is  to  replace  all  the  aluminum 
surfaces  in  the  flow  loop  with  another  material  such  as  stainless  steel  that  will  not 
corrode. 

Another  improvement  to  flow  quality  would  be  the  installation  of  a  deaeration 
system.  At  present,  a  delay  of  two  days  is  required  between  filling  the  channel  and 
making  measurements  in  order  for  the  dissolved  air  to  come  out  of  solution.  Some 
of  this  air  becomes  trapped  in  the  lower  parts  of  the  return  leg,  where  it  is  very 
difflcult  to  remove.  A  deaeraticMi  system  would  reduce  greatly  this  delay  and  would 
also  make  clearer  video  images  possible  for  quantitative  flow  visualization  using 
hydrogen  bubble  timelines  (see  Bruneau,  1992).  It  would  also  make  possible  the  use 
of  hot-film  aimemometry,  which  is  more  accurate  for  making  turbulence 
measurements  than  LDV.  At  this  time  an  effort  to  install  a  deaeration  system  is  being 
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The  quality  of  the  flow  itself  may  be  improved  by  the  final  installation  of  the 
turbulence  control  screens.  Reduction  of  the  particulates  in  the  flow  is  absolutely 
necessary,  however,  before  installation  of  the  turbulence  control  screens  is  possible. 
These  screens  should  reduce  the  turbulrace  level  in  the  freestream,  although  with  the 
present  instrumentation,  the  improvement  may  not  be  measurable. 


CHAPTER  4 


EXPERIMENTAL  RESULTS  AND  DISCUSSION 

4.1  EXPERIMENTAL  CONDITIONS  AND  DATA  PROCESSING 

4.1.1  Experimental  Plan 

To  experimentally  investigate  the  structure  and  developmoit  of  a  turbulmt 
spot,  following  the  objectives  set  forth  in  Chapter  1,  measurements  of  the  unsteady 
velocity  field  both  inside  and  outside  of  the  boundary  layer  were  made  during  the 
passage  of  isolated  turbulent  spots  and  an  interacting  turbulent  spot  pair.  The 
freestream  velocity  during  all  but  one  of  these  experiments  was  U =0.406  m/s.  An 
additional  turbuloit  spot  with  U=0.19  m/s  was  also  investigated. 

Unsteady  streamwise  velocity  profiles  were  measured  in  order  to  calculate  the 
local  fluctuations  in  the  mass  flux  deficit  due  to  the  turbulent  spot.  Five  streamwise 
stations  w^  used  and  were  located  at  20  cm  intmvals,  as  shewn  in  Figure  4.1.  For 
the  isolated  spot  cases  at  U=0.41  m/s,  off-centerline  measurements  were  made  at 
three  spanwise  locations  for  each  streamwise  station.  Unsteady  velocity  profiles  were 
measured  only  on  the  centerline  for  the  interacting  spot  cases.  An  additional  profile 
was  measured  at  x= 1.7S  m  on  the  centerline  for  U  =0. 19  m/s.  The  z-coordinate  for 
the  measuremoits  lay  along  rays  originating  from  the  genoation  point  at  angles  of 
2.5®,  5®,  7.5®  from  the  centerline.  In  addition,  the  component  of  the  velocity 
normal  to  the  wall,  Vp,  was  measured  in  the  potential  flow  above  the  boundary  laym*. 
These  measurements  were  made  at  the  same  (x,z)  locations  as  the  u-component 
measurements,  at  a  height  of  3.67^1,.  An  extra  off-centerline  station  was  added  for 
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Figure  4.1  Top  view  of  locations  for  unsteady  velocity  measurements.  Off- 
centerline  points  on  rays  through  gmeration  point  at  angles  shown,  measured  from 
the  centerline. 


the  v-component  measurements  along  the  ray  through  the  generation  point  at  an  angle 
of  10**  to  the  centerline.  Both  the  isolated  and  interacting  spot  measurements  used 
these  locations  for  measuring  v^.  Hie  locations  and  type  of  measuremoit  for  all 
experiments  are  listed  in  Tables  4.1  -  4.3. 


Table  4.1  List  of  isolated  spot  u-profile  experiments. 


Qfpe  of 

expmmeot 

type  of 

measuieineot 

isolated  q>ot 

u  profile 

isolated  q;x>t 

u  profile 

isolated  qmt 

u  profile 

isolated  qiot 

u  profile 

isolated  spot 

u  profile 

z 

U 

cm 

cm/s 

0,  0.87,  1.74, 
2.63 


0,  1.76,  3.49, 
5.27 


0,  2.62,  5.23, 
7.83 


0,  3.49,  6.97, 
10.53 


0,  4.37,  8.72, 
13.17 


isolated  qiot 


u  profile 


175 


19.1 


Table  4^  List  of  isolated  spot  direct  normal  velocity  measurements. 


type  of 

experiment 

type  of 

measumnent 

X 

m 

z 

cm 

isolated  spU 

direct  v. 

73 

0,  0.87,  1.74, 

2.63,  3.47 

isolated  qiot 

direct  v. 

93 

0,  1.76,  3.49, 

5.27,  6.95 

isolated  spot 

direct 

113 

0,  2.62,  5.23, 

7.83,  10.41 

isolated  spot 

direct 

133 

0,  3.49,  6.97, 

10.53,  13.89 

isolated  qxM 

direct  Vp 

153 

0,  4.37,  8.72, 

13.17,  17.37 

73 


4.1.2  Data  Processing 

The  event  timing  for  the  ensemble  averaging  was  triggered  on  the  spot 
generator  disturbance  pulse.  The  averaged  unsteady  velocity  time  series  formed  the 
basis  for  all  subsequent  calculations.  The  unsteady  streamwise  mass  flux  deficit  for 
a  given  (x,z)  location  was  calculated  by  integrating  the  ensemble-averaged  velocity 
time  series: 


MFD(x,z,t)  -  f /j  _  u(x,y,z,t  )v  . 
PU  I'  U  '  ^ 


(4.1) 


using  the  trapezoidal  rule.  The  velocity  normal  to  the  plate,  v^,  is  related  to 
MFD/pU  according  to: 

v,(x,z,t)  _  a  MFD(x,z,t)  _  3  f  w(x,y,z,t)^  (4^) 

u  ax  pu  az|  u  ^ 


On  the  centerline  of  a  spot  dw/dz  =  0  by  symmetry  (also  shown  by  Wygnansld,  et 
al. ,  1976),  so  Vg  there  is  simply  the  streamwise  derivative  of  MFD/pU.  In  these 
calculations  the  spatial  derivative  in  Equation  4.2  has  been  iq>proximated  by  a 
temporal  derivative  according  to  the  Taylor  hypothesis: 

v.(x,z,t)  -  (4J) 
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The  derivative  was  calculated  using  a  4-point  finite-difference  scheme.  It  should  be 
pointed  out  that  another  choice  for  the  vdocity  in  Taylor’s  hypothesis  is  to  use  the 
convection  speed  of  the  turbulent  spot,  Ue,  which  is  later  shown  to  be  76%  of  the 
freestream  velocity.  Using  U  rather  than  Ue  results  in  a  value  for  v.  which  is  76% 
of  the  value  obtained  if  was  used.  However,  since  the  interest  is  in  the  growth  rate 
and  the  arrival  times  of  the  zeroes  of  v.,  the  absolute  magnitude  is  not  an  issue. 


4.2  SINGLE  SPOT  MEASUREMENTS 
4.2.1  Spot  Growth 

In  order  to  quantify  the  growth  of  the  isolated  turbuloit  spot,  the  edges  of  the 
spot  were  determined  based  on  a  perturbation  velodty  critoia.  The  perturbation 
velodty,  Up,  is  defined  as: 


Up(x,y,z,t) 


u(x,y,z,t)  -  u, 
U 


(4.4) 


The  velodty  u,  corresponds  to  the  velodty  in  the  undisturbed  laminar  boundary  layer 
which  was  determined  by  averaging  the  values  of  the  first  100  time  stq)s  of  the 
velocity  trace,  before  the  arrival  of  the  turbulent  spot.  The  perturbation  velocity  time 
traces  were  thm  low-pass  Altered  by  averaging  Up  at  each  time  stq)  with  the  valu^ 
at  the  neighboring  ±  IS  time  stq)s.  This  averaging  window  correqxmds  to  a 
convolution  of  the  velodty  time  series  with  a  square  pulse  of  unit  hdght  and  of 
0. 124  s  width.  The  velodty  perturbation  information  was  used  solely  to  determine 
the  gross  scaling  of  the  spot  growth  for  comparison  with  other  researchers.  The  edges 
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of  the  spot  were  located  based  on  both  a  level  and  a  slope  criteria.  By  following  the 
Up  time  series  in  the  positive  time  dimction,  the  leading  edge  of  the  spot  was  detected 
by  locating  the  time  step  at  which: 

|uj>0.02  and  |  ^  |  >  3.0 

dt  dt 

where  (dUp/dt).^  is  the  maximum  slope  of  Up(t)  in  the  initial,  laminar  portion  of  the 
time  series.  The  trailing  edge  was  found  by  the  same  method,  but  by  following  the 
time  series  from  the  last  point  and  moving  backwards  in  time.  The  spot  envelope  is 
formed  by  the  locus  of  these  points. 

By  plotting  the  position  of  the  measuremrat  station  vs.  the  time  of  arrival  of 
a  given  portion  of  the  turbuloit  spot  envelope  at  a  streamwise  station,  a  curve  is 
constructed  relating  the  feature’s  position  to  the  time  of  its  arrival.  Such  a  plot  is 
referred  to  as  a  celerity  diagram.  The  slope  dx/dt  of  these  curves  gives  the  celerity 
(convection  velocity)  of  the  feature  of  interest.  The  celerities  for  the  leading  and 
trailing  edges  of  the  turbulent  spot  were  determined  at  the  points  nearest  the  wall. 
The  virtual  origin  for  the  spot  growth  was  also  found  from  this  information.  By 
rescaling  the  spot  envelope  into  similarity  coordinates  (Van  Atta  and  Helland,  1980, 
and  Itsweire  and  Van  Atta,  1983): 

,  -  X,  S,  .  0.37(JL)’(i  -  0.3933m)\  {  -  „ ^ ,  (4.«) 

^0*  U^(t  -  tj,) 

it  can  be  see  in  Figure  4.2  that  for  the  last  three  stations  (x  =  1.13  m,  1.33  m,  and 
1.53  m)  the  envelope  collapses  onto  one  curve,  while  the  envelopes  for  the  tirst  two 
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stations  (x  =  0.73  m,  0.93  m)  do  not.  This  result  is  in  agreenumt  with  the  findings 
of  Schubauer  and  Klebanoff  (1956)  and  Wygnanski  et  a/.  (1976)  who  showed  that  the 
spot  settles  into  a  self-similar  structure  some  distance  downstream  of  the  spot 
generator.  This  result  also  indicates  that  the  present  data  include  portions  of  the  initial 
growth  of  the  spot  upstream  of  where  similarity  scaling  ^lies. 

4.2.2  Unsteady  Mass  Flux  Deflcit,  Normal  Velocity 

Figure  4.3  shows  a  typical  trace  of  AMFD/pU  and  v^/U  versus  time. 
Superimposed  on  these  traces  is  the  2%  velocity  perturbation  contour  (the  locus  of 
points  for  which  the  perturbation  velocity  Up,  described  in  Section  4.2.1,  is  ±2%) 
wiiich  shows  the  outline  of  the  turbulait  spot  relative  to  the  features  of  AMFD/pU  and 
Vg.  At  the  time  of  arrival  of  the  spot  the  value  of  AMFD/pU  increases  rapidly  to  its 
highest  peak.  It  then  reverses  direction  and  Ms  below  zero  to  a  negative  peak.  After 
this  peak,  the  curve  slowly  s^ro^hes  the  initial  state.  A  comparison  of  the 
AMFD/pU  curve  and  the  spot  envelope  helps  to  explain  this  behavior.  The  maximum 
peak  corresponds  to  the  arrival  of  the  maximum  height  of  the  ^t,  while  the 
minimum  peak  corresponds  to  the  trailing  edge.  The  slow  return  to  the  initial  state 
represents  the  "calmed  region",  first  idmtified  by  Schubauer  and  Kld)anoff  (1956). 
This  form  is  concordant  with  results  shown  by  Wygnanski  et  p/.(1976)  and  Antonia 
et  a/.  (1981).  Measurements  taken  at  off-coitm’line  locations  show  similar  trends,  but 
with  a  differ^t  phase  and  amplitude  (Figure  4.4).  At  off-centerline  locations  6*  is  not 
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Figure  4.3  Mass  flux  deficit  and  normal  velocity  fluctuaticms  at  x= 1. 13  m,  z»0, 
U«,=40.6  cm/s.  Spot  edge  is  also  shown,  denoted  by  2%  velocity  disturbance 
contour. 


MFD/pU  (cm),  v,/U 
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well  represented  by  MFD/pU  because  there  is  an  appreciable  spanwise  component  of 
velocity  of  up  to  0.08  U,  according  to  Wygnanski  et  a/.  (1976). 

One  of  the  unr^lved  issues  concerning  the  character  of  the  sound  radiation 
from  a  transitional  flow  is  the  question  of  whetha*  the  flow  is  a  monopole  source  or 
not.  If  it  is  the  net  mass  flux  normal  to  the  wall  caused  by  a  tuibulrat  spot  should  be 
nonzero  for  some  instant.  However,  it  is  easy  to  argue  that  the  normal  mass  flux 
fluctuation  due  to  the  spot  must  always  be  zero.  If  we  consider  two  y-z  planes, 
infinite  in  the  z-diiection,  far  enough  upstream  and  downstream  of  a  turbuloit  spot 
that  the  flowfield  fluctuation  due  to  the  spot  is  not  felt,  it  is  clear  that  the  velocity 
profiles  at  these  stations  are  those  of  the  undisturbed  laminar  boundary  layer.  The 
amount  of  mass  flux  normal  to  the  wall  between  these  stations,  which  includes  the 
effect  of  the  turbulent  spot,  should  be  equal  to  the  amount  lost  to  the  streamwise  flow 
between  these  stations  due  to  the  growth  of  the  laminar  boundary  layer.  In  other 
words,  the  net  normal  mass  flux  is  the  same  as  if  the  turbulent  spot  were  not  present, 
or,  that  the  net  normal  mass  flux  fluctuation  due  to  the  presence  of  a  turbulent  spot 
is  zero.  Again,  this  demonstrates  conclusively  that  the  monopole  strmgth  of  a 
turbulent  spot  is  identically  zero.  An  equivalent  statement  is  that  the  integral  of  the 
mass  flux  deficit  over  the  spot  area  at  any  instant  is  zero. 

Knowing  this  property  of  the  velocity  field,  a  test  of  the  measurement  of  the 
mass  flux  deficit  fluctuation  was  carried  by  integrating  AMFD/pU(x,z,t)  over  x  and 
z  a.  1  single  instant.  Because  the  mass  flux  deficit  was  measured  at  four  (x,z) 
locations  as  a  function  of  time,  it  was  necessary  to  make  use  of  a  similarity 
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transformation  of  the  data  in  order  to  obtain  the  mass  flux  deficit  distribution  over 
space  for  one  instant.  This  was  accomplished  using  the  similarity  transformation 
given  by  Equations  4.6.  Data  from  the  four  velocity  profiles  taken  at  x = 1 .33  m  were 
used  to  calculate  the  unsteady  mass  flux  deficit.  The  integral  over  the  x-z  plane  was 
transformed  to  coordinates  and  the  integral  calculated  using  the  Trapezoid  Rule 
in  two  dimensions.  The  result  of  the  integration  was  0.8%  of  the  incoming  mass  flux, 
indicating  that  the  mass  flux  deficit  calculations  based  on  the  velocity  measuremoits 
are  accurate  to  that  degree. 

Figure  4.S  shows  die  mass  flux  deficit,  normal  velocity  fluctuations,  and  2% 
perturbation  velocity  contour  for  the  profile  taken  at  x=1.75  m,  U=0.19  m/s.  Note 
that  the  amplitudes  of  the  peaks  are  larger  than  those  for  the  higher  speed  case, 
consistent  with  the  thicker  boundary  layer  preset  in  a  lower  speed  flow.  In  addition, 
the  total  time  for  a  spot  passage  is  much  larger.  The  overall  sh^  of  both  the  mass 
flux  deficit  and  the  normal  velocity  are  similar  to  the  higher  speed  case,  but  the  peak 
values  of  \J\J  are  larger  and  the  decay  of  the  mass  flux  deficit  from  the  positive 
maximum  to  the  negative  minimum  peak  is  much  sharper.  In  addition,  the  small 
negative  peak  in  \J\J  pnaent  in  the  higher  speed  case  is  barely  visible  in  Figure  4.5. 
More  discussion  of  the  difierences  between  these  two  cases  is  given  in  Section  4.2.3. 

Further  insight  into  the  behavior  of  the  unsteady  AMFD/pU  may  be  gained  by 
inspecting  Figure  4.6,  which  shows  the  time  series  for  AMFD/pU  and  representative 
time  series  for  the  u-velocity  in  the  boundary  layer.  Deep  in  the  boundary  layer  the 
velocity  undergoes  a  sudden  acceleration  upon  the  arrival  of  the  spot  leading  edge. 
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Figure  4.5  Mass  flux  deficit  and  normal  velocity  (vj  fluctuations  for  X-1.7S  m, 
z=0  m,  U=19  cm/s.  Also  shown  is  the  spot  edge,  denoted  by  the  2%  velocity 
perturbation  contour. 
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Figure  4.6  Comparison  of  velocity  pertuibation  contours  with  mass  flux  deficit 
time  series,  x=1.13  m. 
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At  greater  distances  above  the  wall  this  acceleration  becomes  less  pronounced  and  is 
limited  to  the  arrival  of  the  spot  leading  edge,  after  which  the  flow  decelerates.  The 
initial  velocity  acceleration  in  the  region  corresponding  to  the  undisturbed  laminar 
boundary  layer  arrives  at  earlier  times  as  the  distance  from  the  wall  increases.  The 
mass  flux  deficit  increases  from  the  undisturbed  value,  but  only  slightly,  in  this  region 
of  the  spot  which  has  been  called  the  "overhang”  r^ion  by  Wygnanski  et  o/.(1976). 
The  deceleration  of  the  flow  following  the  overhang  is  related  to  the  increase  in  the 
boundary  layer  height  in  the  central  part  of  the  spot.  Here  the  spot  actually  protrudes 
above  the  height  of  the  surrounding  laminar  boundary  layer. 

The  deceleration  of  the  flow  outside  the  boundary  layer  corresponds  exactly 
to  the  sharp,  positive  peak  in  the  mass  flux  deficit.  This  maximum  in  turn 
corresponds  not  only  to  the  minimum  in  the  u-velocity  peaks  but  also  to  the  maximum 
spot  height.  After  this  maximum,  the  velocity  increases  again.  The  flow  which  was 
initially  outside  the  laminar  boundary  layer  quickly  accelerates  to  its  undisturbed 
value.  The  flow  in  the  boundary  layer  continues  to  accelerate  to  a  maximum  which 
corresponds  to  the  trailing  edge  of  the  turbulent  part  of  the  spot  structure.  This  was 
also  (^served  by  Wygnanski  et  fi/.(1976).  The  mass  flux  deficit  reaches  its  minimum 
at  this  point.  Following  the  arrival  of  the  trailing  edge  of  the  turbuloit  r^on,  the 
flow  is  again  laminar  and  it  slowly  decays  to  its  undisturbed  level. 

The  normal  velocity  on  the  spot  centerline,  calculated  from  the  AMFD/pU,  is 
also  shown  for  x  ==  1 . 13  m  in  Figure  4.3.  First,  note  that  the  normal  velocity  initially 
undergoes  a  r!q)id,  negative  fluctuation,  foUowed  by  a  positive  peak  which  is  at  first 
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sharp  and  then  more  gradual.  The  positive  peak  is  in  turn  followed  by  a  shallow 
negative  peak  where  the  velocity  slowly  returns  to  its  initial  value.  The  two  larger 
leading  edge  peaks  are  seen  to  be  associated  with  the  positive  peak  in  MFD/pU,  while 
the  normal  velocity  fluctuations  near  the  spot  trailing  edge  are  much  smaller. 

The  spanwise  variation  in  v,  is  shown  in  Figure  4.7.  Note  how  the  negative 
peak  at  the  spot  leading  edge  occurs  at  later  times  as  the  distance  from  the  craterline 
increases.  For  comparison  the  fluctuating  velocity  component  normal  to  the  wall,  Vp, 
is  also  shown  in  Figure  4.7.  It  was  measured  directly  at  a  height  of  3.675i.  above  the 
wall.  It  is  important  to  note  that  Vp  is  defined  as  the  velocity  normal  to  the  wall  minus 
the  steady  component  of  normal  velocity  due  to  the  growth  of  die  boundary  layer. 
The  overall  shape  of  the  time  traces  for  v.  and  Vp  are  similar.  However,  the  phase 
of  the  curves  is  very  different,  especially  when  sera  from  a  comparison  of  the 
spanwise  distributions.  The  phase  of  Vp  is  preserved  across  the  span  of  the  spot,  while 
for  Vg  it  is  not.  Note  that  for  Vp  the  negative  leading  edge  peak  disappears  as  the 
distance  from  the  centerline  increases. 

Contour  plots  of  Vp  at  x=0.73  m,  x=1.13  m,  and  x=1.53  m  are  shown  in 
Figures  4.8-4.10.  The  development  of  the  potential  flow  disturbance  can  be  sera 
clearly.  Not  only  do  the  peaks  in  velocity  grow  in  amplitude  (the  negative  leading 
edge  peak  in  particular),  but  the  area  covered  by  each  peak  also  increases  with 
downstream  distance.  This  spatio-temporal  bdiavior  is  important  for  the  prediction 


of  the  radiated  sound. 
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Figure  4.7  Comparison  of  normal  velocity  deduced  from  mass  flux  deficit 
fluctuations  to  those  obtained  by  direct  measurement,  x=1.13  m. 
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Figure  4.8  Contour  plot  of  perturbation  velocity  normal  to  the  plate  measured  at 
a  height  of  3.675^  for  x=  0.73  m,  U=40.6  cm/s.  Dotted  lines  denote  n^adve 
velocity,  solid  lines  positive. 
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The  discrepancy  in  the  behavior  between  v„  and  Vp  can  be  explained  by 
reference  to  Equation  2.3,  where  the  freestream  velocity  potoitial  fluctuations  are 
expressed  in  terms  of  the  integral  of  the  normal  velocity  fluctuations  due  to  unsteady 
displacement  thickness  fluctuations.  The  normal  derivative  of  this  equation  yields  the 
relationship  betwe^  Vp  and  Vp.  For  this  situation  where  the  observer  is  in  close 
proximity  to  the  source,  the  retarded  time  may  be  neglected  because  r/c-K),  resulting 
in  the  equation  given  by  Van  Atta  et  a/.  (1982): 


Vp(x,t) 


p  r  rv.(x,z,t)(y'  -  y) 

2rll 


(4.7) 


They  showed  experimentally  that  the  value  of  Vp  on  the  centerline  was  dq)endent  on 
the  entire  boundary  layer  flowfield  through  Equation  4.7.  From  this  equation  it  is 
apparent  that  for  a  given  instant  the  p^tuibation  velocity  in  the  freestream  is 
influoiced  by  the  entire  boundary  layer  flowfield  at  that  instant,  weighted  by  the  third 
power  of  the  distance  from  the  "source"  point  in  the  boundary  layer  and  the 
measurement  point  in  the  freestream.  It  is  also  necessary  to  recall  that  the  method 
used  here  for  calculating  Vp  does  not  consider  the  spanwise  component  of  velocity  (see 
Equation  4.3). 

The  streamwise  variation  of  the  unsteady  mass  flux  deficit,  measured  at  five 
X  locations,  is  shown  in  Figure  4.12.  The  peak  values  clearly  increase  with  x,  while 
the  peaks  themselves  ^)read  apart  and  become  less  steq>.  The  streamwise  variation 
in  the  three  parameters  of  interest,  the  magnitude  of  the  mass  flux  deficit  fluctuations. 
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Figure  4.11  Streamwise  variaticm  of  unsteady  mass  flux  deficit,  z=0,  Ug,=40.6 
cm/s.  Rise  time  tj  is  defined  as  the  difference  in  times  of  arrival  of  leading  edge 
and  first  positive  peak. 
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the  spot  convection  velocity,  and  the  rise  time  may  all  be  deduced  from  the 
information  in  Figure  4. 10.  Of  considerable  interest  in  analyzing  the  acoustics  is  the 
streamwise  growth  of  the  magnitude  of  the  fluctuations  of  AMFD/pU  as  the  spot 
convects  downstream.  Figure  4.11  shows  both  the  height  of  the  positive  peak  and  the 
difference  in  heights  between  the  maximum  and  minimum  of  AMFD/pU  vs. 
streamwise  location. 

A  celerity  diagram  of  the  principle  features  of  the  local  AMFD/pU  curve  is 
shown  in  Figure  4.13.  The  location  of  each  measurement  is  plotted  against  the  time 
of  arrival  of  each  of  these  features.  The  time  of  arrival  of  the  leading  edge  was 
determined  in  a  manner  similar  to  that  previously  described  which  was  used  to  locate 
the  arrival  of  the  edges  of  the  turbulent  spot,  as  described  above.  The  time  of  arrival 
of  the  positive  and  negative  peaks  was  found  by  locating  the  zeroes  of  the  time 
derivative  of  the  AMFD/pU  curve.  A  straight-line  curve  was  fit  to  the  data  from  the 
four  downstream  stations.  These  curves  are  also  shown  in  Figure  4.13.  Notice  that 
these  curves  do  not  intersect  at  the  spot  generation  location.  This  result  is  similar  to 
the  findings  of  other  investigators  and  can  be  explained  by  the  initial  nonlinear  growth 
of  the  turbulent  spot.  The  convection  speed  for  the  given  feature  of  interest  is  given 
by  the  slope  of  the  straight-line  curve.  Hie  spot  arrival  moves  with  a  convection 
speed  of  U,  the  maximum  peak  with  speed  0.76  U,  and  the  minimum  peak  with  speed 
0.52  U.  The  last  two  convection  speeds  correspond  well  with  the  overall  and  trailing 
edge  convection  speeds  of  the  isolated  spot  observed  by  Cantwell  er  p/.(1978). 
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Figure  4.13  Celerity  diagram  for  MFD/pU  shown  in  Fig.  7.  Convection  speeds 
are  U  for  leading  edge,  0.76  U  for  (+)  peak,  and  0.52  U  for  (-)  peak. 
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The  streamvase  variation  in  the  peak  amplitudes  of  the  two  large  normal 
velocity  peaks  are  shown  in  Figure  4.14.  Both  the  positive  and  negative  velocity 
peaks  increase  with  x,  levelling  off  at  x  •>  1. 13  m  from  the  plate  leading  edge.  The 
most  intense  normal  mass  flow  fluctuation  occurs  downstream  of  this  location.  The 
slopes  of  these  curves  are  greatest  in  die  most  upstream  portion  of  the  intermittrat 
region,  indicating  that  the  rate  of  change  of  the  normal  mass  flux  is  greatest  there. 


4.2.3  Rise  Time  Estimates 

The  AMFD/pU  time  traces  clearly  show  that  the  rise  and  fall  times  which 
characterize  the  large-scale  fluctuation  in  displacement  thickness  are  very  different 
from  one  another.  These  two  parameters,  which  we  name  t*  (rise)  and  r  (fall),  are 
deduced  from  the  celerity  diagrams.  The  rise  time  for  a  given  (x,  z)  location  is 
considered  to  be  equal  to  the  time  at  which  the  maximum  peak  arrives  minus  the  time 
at  which  the  spot  arrives  (Figure  4. 10).  Similarly,  the  fall  time  is  defined  as  the 
difference  in  the  arrival  time  of  the  maximum  and  minimum  peaks.  The 
nondimensional  rise  time,  Uet^/Ax,  identified  by  Lauchle  (1981)  as  controlling  the 
high-frequency  content  of  the  radiated  sound  spectrum,  may  be  estimated  using  these 
parameters.  However,  in  order  to  estimate  the  non-dimensional  rise  time  a  transition 
zone  loigth  is  needed.  An  estimate  of  the  transition  zone  length  Ax  was  obtained 
using  the  correlation  (Josserand  and  Lauchle,  1990): 

Re^  -  17ReJ‘  (4.8) 

with  Xo  =  .53  m  (^t  generator  location).  This  correlation  predicts  Ax  •  77  cm  for 
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U=40.6  cm/s,  and  Ax  •  90  cm  for  U= 19  cm/s.  A  corresponding  normalized 
transition  zone  streamwise  coordinate,  f,  may  be  defined  by: 

^  (X  -  X,)  (4 

*  Ax 

Figure  4.14  shows  how  UeVAx  increases  with  f.  The  ranges  of  the  normalized  rise 
and  fall  times  for  the  40.6  cm/s  spot  are: 

0.06  <  <  0.30  0.13  <  ^  <  0.74  W-IO) 

Ax  Ax 

Figure  4.15  shows  a  clear  monotonic  increase  in  the  magnitude  for  both  the  rise  and 
fall  times.  From  these  results  it  can  be  seen  that,  for  an  isolated  turbulent  spot,  the 
rise  time  of  the  intermittent  fluctuation  is  always  less  than  the  time  for  the  spot  to 
traverse  the  intermittent  zone.  This  is  especially  true  for  the  upstream  portions  of  the 
intermittent  zone.  The  increase  in  the  rise/fall  times  is  due  to  the  growth  of  the 
streamwise  extent  of  the  qrot  and  seems  to  indicate  that  the  mass  flux  normal  to  the 
plate  decreases  with  x.  As  shown  in  Figure  4. 10,  the  magnitude  of  the  displacement 
thickness  fluctuation  increases  with  x,  which  implies  a  more  intense  local  fluctuation 
in  V..  The  balance  between  these  two  effects  explains  the  levelling  off  of  the  peak 
normal  velocity  shown  in  Figure  4.13. 

Pmaud  (1989)  compared  an  estimate  of  the  wall  pressure  qrectrum  in  the 
laminar  portion  of  the  boundary  layer  upstream  of  the  transititm  with  measured 
wall  pressure  spectra.  The  only  pressure  disturbances  presumed  to  be  measured  in 
the  laminar  portion  of  the  boundary  layer  are  those  which  propagate  as  sound  from 
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other  regions  of  the  flow.  The  experiments  were  carried  out  on  a  body  of  revolution 
in  air.  The  spectral  estimates  were  made  using  the  ai^roach  of  Lauchle  (1981)  for 
T=0.05  and  T=0.1S.  These  estimates  were  found  to  bracket  the  measured  spectrum 
well,  except  for  some  low-ftequency  peaks  which  were  attributed  to  ToUmien- 
Schichting  waves.  The  range  of  T  from  the  study  of  Perraud  (1989)  agrees  well  with 
that  found  in  our  measurements  of  the  rise  time.  This  correspondence  suggests  that 
the  flow  measured  here  should  exhibit  similar  pressure  behavior. 

Also  shown  in  Figure  4. 14  are  the  rise  and  fall  times  for  the  mass  flux  deficit 
of  a  turbulent  spot  at  11=0.19  m/s,  x=1.75  m  from  the  leading  edge.  Note  that  T^ 
seems  to  follow  the  traid  for  the  higho*  vdocity  case,  while  T  does  not.  This  result 
may  indicate  that  T'^  is  insensitive  to  the  Reynolds  number.  Wygnanski  et  a/.  (1982) 
showed  that  the  streamwise  growth  rate  of  a  turbulent  spot  was  sensitive  to  the 
Reynolds  number  based  on  the  undisturbed  flow  displacement  thickness  at  the 
generator  location,  (Re^.)^.  For  the  two  cases  given  here,  (Re«.)^  was  equal  to  S37 
and  797  for  U  =  0.19  m/s  and  0.41  m/s,  respectively.  The  ratio  of  the  trailing  edge 
convection  speed  to  the  freestream  speed  was  found  to  decrease  with  an  increase  in 
(Rc<*)tai-  However,  the  leading  edge  and  spot  maximum  heights  w^  found  to  be 
insensitive  to  Reynolds  number.  The  spot  leading  edge  and  maximum  hdght 
correspond  to  the  leading  edge  and  positive  maximum  of  the  mass  flux  defidt,  while 
the  spot  trailing  edge  corresponds  to  the  negative  mass  flux  defidt  peak,  as  shown  in 
Figure  4.S.  Since  the  rise  and  M  tim^  are  related  to  the  convection  speeds  of  these 
features  of  the  mass  flux  defidt,  it  follows  that  the  nondimensional  rise  time  should 
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be  insensitive  to  (Rea.)gaif  while  the  nondimensional  fall  time  should  decrease  with 
(Rea*)(ea>  which  is  agreement  with  the  observations  presented  here.  Further 
investigation  should  more  strongly  substantiate  this  assertion. 

4.3  INTERACTING  SPOT  MEASUREMENTS 
4.3.1  Unsteady  Mass  Flux  Deficit,  Normal  Velocity 

Figure  4. 16  shows  the  streamwise  development  in  the  mass  flux  deficit  on  the 
centerline  for  the  int^acting  spot  case.  The  bdhavior  of  each  spot  at  the  upstream 
stations  is  much  the  same  as  for  two  indqiendrat  spots,  excq)t  that  the  mass  flux 
deficit  positive  peak  for  the  second  spot  has  a  reduced  amplitude.  As  the  second  spot 
catches  the  first,  the  growth  of  the  mass  flux  deficit  positive  peak  is  attenuated  relative 
to  the  peak  for  the  first  spot.  The  second  spot  leading  edge  and  the  first  spot  trailing 
edge  are  just  meeting  at  x=1.33  m,  and  have  merged  at  x=l.S3  m. 

Figures  4.17  and  4.18  ^ow  the  relation  between  the  unsteady  velocity  field 
and  the  mass  flux  deficit  for  x=0.93  m  and  x=l.S3  m.  Recalling  the  discussion 
concerning  Figure  4.S,  these  comparisons  show  how  the  presence  of  the  first  spot 
affects  the  developmoit  of  the  velocity  field,  and  hence  the  mass  flux  deficit,  of  the 
second  spot.  For  the  case  x=0.93  m,  the  two  spots  are  seoi  to  be  distinct,  although 
the  leading  edge  of  spot  2  has  oicroached  on  the  calming  region  of  spot  1.  The 
AMFD/pU  positive  peak  of  ^t  2  is  attenuated  compared  to  that  of  spot  1.  The 
corresponding  attenuation  of  the  velocity  minima  in  the  uiq)er  r^ons  of  the  qx>t  may 
be  seen  in  Figure  4. 16.  Another  way  of  stating  this  result  is  that  the  second  spot  has 
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Figure  4.17  Comparison  of  the  velocity  Held  inside  the  interacting  turbulent  ^ts 
and  the  correspcmding  mass  flux  deficit  for  x=0.93  m.  The  isolated  spot  mass  flux 
deficit  is  shown  for  lefnence. 
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Figure  4.18  Compariscm  of  velocity  traces  with  mass  flux  for  interacting  spots  at 
x=  1.53  m.  The  leading  edge  of  spot  2  has  m^ed  with  the  trailing  edge  of  ^t 
1,  as  indicated  by  the  dissy)pearance  of  the  overhang  region  in  spot  2. 
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a  lower  maximum  height  than  the  first.  This  situation  is  even  more  pronounced  in 
Figure  4.17,  which  shows  the  interaction  for  x=l.S3  m.  Here  it  is  clear  that  the 
initial  acceleration  of  the  flow  at  the  spot  leading  edge,  which  characterizes  the  spot 
overhang  region,  has  vanished  in  spot  2  due  to  its  moging  with  spot  1.  Although  the 
spots  have  merged  at  this  location,  the  two  spots  still  have  distinct  structure:  the 
velocity  minima  and  maxima  are  still  clearly  identifiable,  as  are  the  corre^nding 
features  of  AMFD/pU,  albeit  with  some  modification  in  amplitude.  Whether  or  not 
the  identities  of  the  spots  would  become  blurred  further  downstream  is  not  clear  at  this 
point  and  motivates  further  study. 

A  relevant  question  is  how  the  spot  convection  velocity  is  modified  by  the 
interaction.  To  answer  this  question  celerity  diagrams  for  the  interacting  spots  were 
constructed  using  the  method  described  in  Section  4.2.2.  These  are  shown  in  Figure 
4. 19.  The  spot  convection  speed  is  given  by  the  slope  of  the  curve  corresponding  to 
the  positive  AMFD/pU  peak.  It  is  clear  that  the  convection  q)eed  has  not  been  altoed 
significantly  by  the  interaction. 


4.3.2  Rise  Time  Estimates 

Rise  time  estimates  for  the  interacting  spot  case  were  made  as  described  above 
in  Section  4.1.2.  Figure  4.20  shows  the  results.  The  same  trends  and  ordo^  of 
magnitude  that  held  for  the  isolated  spot  case  are  clearly  seen.  These  results  are  not 
surprising  given  that  the  rise  and  fall  times  are  defined  as  the  difference  between  the 
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Figure  4.20  Celerity  diagrams  for  interacting  ^t  cases.  Spot  1  is  the  first 
generated.  2  doiotes  second  spot  and  i  draotes  isolated  spot  case.  Arrival  time 
defined  relative  to  generation  time. 
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arrival  times  of  the  various  features  of  the  spot  and  these  have  not  been  seriously 
modified  by  the  interaction  between  the  two  turbulent  spots. 


CHAPTERS 


ESTIMATION  OF  SOUND  RADIATION 

The  measurem^its  described  in  Chapter  4  may  be  used  to  estimate  the  direct 
acoustic  emissions  from  a  single  turbulent  spot  and  a  weakly  interacting  spot  pair. 
The  results  of  these  cases  may  also  be  used  to  estimate  the  sound  radiation  from  a 
naturally  occurring  transition  region  after  making  some  simplifying  assumptions.  The 
estimation  of  the  sound  radiation  is  carried  out  in  the  spirit  of  the  Liepmann  analogy, 
where  it  is  related  to  the  unsteady  normal  velocity  distribution  on  a  surface  outside  the 
boundary  layer.  By  observing  the  effect  of  a  turbulent  spot  passage  on  the  normal 
velocity  distribution  on  this  surface,  some  conclusions  may  be  drawn  about  the  nature 
of  the  acoustic  source.  These  observations  and  the  experimental  data  may  thoi  be 
used  to  provide  a  quantitative  estimate  of  the  sound  radiation  from  a  turbulent  ^t. 
Finally,  the  assumption  that  the  displacemoit  thickness  is  realistically  rqnesented  by 
an  intermittency-weighted  form  hypothesized  by  Lauchle  (1981)  and  others  is 
evaluated. 

5.1  SOUND  RADIATION  FROM  AN  ISOLATED  TURBULENT  SPOT 

In  considering  the  sound  radiation  from  a  single  turbulent  spot,  two  general 
statements  may  be  made  concerning  the  nature  of  the  acoustic  source  by  inspection  of 
the  data  presented  in  ChsqMer  4.  The  first  is  that  the  source  is  compact,  givm  the  low 
flow  Mach  number.  For  this  reason,  the  retarded  time  differences  x7c  in  the 
argument  of  the  int^rand  of  Equation  2.6  may  be  neglected. 
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The  second  statement  that  can  be  made  about  the  nature  of  the  acoustic  source 
follows  from  the  incompressibility  of  the  source  flow.  Incompressibility  dictates  that 
the  net  instantaneous  mass  flux  through  a  surface  parallel  to  the  wall  must  be  equal 
to  the  amount  of  mass  flux  lost  to  the  incoming  stream  due  to  steady  growth  of  the 
laminar  boundary  layer.  In  other  words,  the  integral  of  Vp  over  the  surface  parallel 
to  the  wall  is  zero  at  all  times.  At  any  x-location  the  distribution  of  unsteady  mass 
flux  normal  to  the  plate  consists  of  three  peaks:  a  sharp  negative  peak  at  the  spot 
leading  edge,  followed  by  a  large  positive  peak  whose  trailing  edge  decays  slowly, 
and,  finally,  a  weak  negative  peak  at  the  spot  trailing  edge.  As  the  spot  develops, 
each  of  these  peaks  spreads  and  the  maximum  velocity  increases.  Therefore,  the  form 
of  the  unsteady  Vp  vs.  time  curve,  as  shown  in  Figures  4.6-4.9,  suggests  a  multipole 
source  mechanism. 

The  unsteady  mass  flux  distribution  through  this  surface  due  to  the  turbulent 
spot  passage  may  be  modeled  as  the  superposition  of  four  simple  time-dq)endent 
sources  as  shown  in  Figure  5.1.  The  leading  edge  consists  of  two  sources  of  equal 
strength  Q,(t)  but  opposite  phase,  sqxarated  by  a  distance  6=Ueti.  The  trailing  edge 
consists  of  two  weaker  sources  of  strength  Q2(t)  similarly  arranged.  Thus,  the  spot 
appears  to  be  the  sum  of  two  unequal  dipoles  of  opposite  phase.  The  leading  edge 
dipole  may  be  further  decomposed  into  two  coincident  dipoles,  one  composed  of 
simple  sources  of  strength  Qt-Q2  and  the  other  of  sources  of  strength  Qj.  This 
decomposition  allows  the  spot  to  be  modeled  as  the  sum  of  a  dipole  of  strragth  cd(Qi- 
Q^/dt  and  a  quadrupole  of  strength  €^dQ2/dt.  Because  the  spot  increases  in  length  as 
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Figure  5.1  Spatial  distribution  of  instantaneous  v.  due  to  a  turbulent  spot  on  a 
surface  outside  the  boundary  layer  and  parallel  to  the  wall,  and  its  equivaloit 
multipole  source  distribution. 
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it  develops,  the  peaks  spread  apart.  Thus,  the  source  strength  will  vary  not  only 
because  dQ/dt  changes  but  also  because  6  increases. 

Having  made  these  qualitative  statemoits  about  the  nature  of  the  acoustic 
source,  the  expmimental  results  presoited  in  Chapter  4  may  be  used  in  a  scaling 
analysis  of  Equation  2.6  to  give  some  quantitative  estimates  of  the  sound  radiation. 
Before  the  scaling  is  carried  out,  the  integral  in  Equation  2.6  must  be  rewritten  in  a 
way  that  accounts  for  the  source  compactness  and  for  the  multipole  nature  of  the 
source  flowfield.  If  this  is  not  done,  a  scaling  analysis  based  on  the  compact  form  of 
Equation  2.6  will  result  in  an  overestimation  of  the  sound  level,  since  this  form  of  the 
integral  thoi  describes  a  monopole  source,  and  it  has  been  established  that  the 
monopole  straigth  is  zero.  In  order  to  account  for  phase  differences  in  the  compact 
source  distribution.  Equation  2.6  may  be  expressed  as  a  multipole  expansion  (see 
Pierce  (1989)): 


P(r,t)  * 
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Ids'  +  1! 


(5.1) 

where  [  ]  denotes  that  the  argument  is  evaluated  at  the  retarded  time  t-r/c  and  jSj  is  a 
directivity  factor.  As  discussed  above,  the  first  term  describing  monopole-type 


radiation  is  zero. 
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Using  the  following  scaling: 


V  — 


-  Ut 

C  % 


2t2 


ua 

C  I 


(5^) 


the  dipole  radiation  pressure  pa,  source  stmigth  D,  and  radiated  power  IPa  scale  as: 
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The  quadrupole  radiation  pg,  source  stroigth  q,  and  radiated  power  !Pq  scale  as: 


Pq  -  ^H^M^Aa  -  /S^pj  q 


-  MD 


(5.4) 


(Pg  ~  M*(P, 

The  relative  powers  of  Mach  number  indicate  that  the  quadrupole  radiatimi  from  a 
turbulent  spot  is  n^ligible  compared  to  the  dipole  radiatitm  at  low  Mach  numbos. 
From  this  scaling  argument  it  can  be  ccmcluded  that  the  sound  radiatitm  from  the 
large-scale  motion  in  a  turbulent  spot  is  dominated  by  the  dipole  associated  with  the 
positive  peak  of  the  mass  flux  deficit  distribution. 

It  ^ould  be  pointed  out  that  using  the  spoi  size  as  a  length  scale  is  equivaloit 
to  assuming  that  Vp  is  correlated  over  the  q[)Ot.  This  assumption  is  justified  by 
reference  to  the  oisemble-avoaged  measurements  of  Vp,  which  show  a  cohoent 
velocity  pattern  over  the  spot.  It  was  previously  shown  by  Van  Atta  et  al.  (1982)  that 
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the  fluctuations  in  the  potential  flow  velocity  are  repeatable  and  nonturfoulent, 
suggesting  that  the  ensemble  average  merely  filters  out  noise  associated  with  the 
measurement,  rather  than  averaging  out  turbulent  fluctuations. 

Using  the  experimental  results  presmted  above,  the  streamwise  variation  of  the 
radiated  sound  level  may  be  deduced.  Figure  5.2  shows  the  variation  of  D  and  p^ 
with  X.  Because  D  is  proportional  to  \l\^  the  dipole  source  strength  per  unit  area  is 
higher  in  the  early  stages  of  the  spot  development.  The  maximum  is  not  at  the 
furthest  upstream  station,  however.  This  effect  is  due  to  the  low  amplitude  of  the 
fluctuation  in  the  mass  flux  deficit  at  the  upstream  stations.  Because  the  source  area 
Amukc  is  proportional  to  dq)«idrace  of  pj  on  t;  cancels  and  the  streamwise 

variation  in  radiated  sound  depends  oitirely  on  the  growth  in  AMFD/pU  with  x.  The 
transition  noise  calculations  of  Lauchle  (1981)  and  Audet  er  a/.  (1989)  did  not  include 
the  streamwise  variation  in  either  T  or  the  peak  value  of  AMFD/pU. 

5.2  INTERMITTENCY- WEIGHTING  OF  DISPLACEMENT  THICKNESS 

The  form  of  the  measured  MFD/pU  differs  somewhat  from  that  assumed  by 
Lauchle  (1981)  as  indicated  by  Figure  5.3.  First,  the  measured  level  within  the  spot 
is  not  constant.  Second,  AMFD/pU  takes  on  negative  values  in  the  portions  of  the 
spot  which  arrive  last  at  the  measurement  location.  This  second  difference  is  the  most 
important,  since  it  reflects  the  principle  of  conservation  of  mass.  Because  the  source 
flow  is  incompressible,  the  net  mass  flux  through  So  at  any  instant  is  identically  zero. 
For  the  assumed  bdiavior,  the  growth  of  the  spot  will  result  in  the  leading  and  trailing 


0.6  0.7  0.8  0.9  1.0  1.1  1.2  1.3  1.4  1.5  1.6 

x(m) 

Figure  5.2  Order  of  magnitude  estimates  of  the  radiated  pressure,  p^,  source  area, 
and  dipole  strength  per  unit  area,  D,  from  an  isolated  turbident  qxrt,  vs. 
streamwise  coordinate  for  U =0.406  m/s. 
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Figure  5.3  Coinparis(Hi  of  measured  MFD/pU.  and  intermittracy-weighted  d*(t) 
assumed  by  Lauchle  (1981).  Assumed  form  modds  only  positive  peak  associated 
with  dominant  dipole  source. 


edges  of  the  bulge  in  5*(x,t)  spreading  aq)art.  This  effect,  coupled  with  the  growth  in 
the  magnitude  in  A5*  with  x,  which  both  Lauchle  (1981,  1989)  and  Audet 
a/.  (1989a,  b)  n^lected,  indicates  that  the  volume  bounded  by  the  solid  wall  and  the 
surface  y  =  d*(x,z,t)  is  growing  slowly  with  time.  This  cannot  occur  if  mass  is  to  be 
conserved  in  the  flow.  If  A5*(x,z,t)  and,  hence,  AMFD/pU  has  a  positive  peak  at  one 
location,  it  must  also  have  a  negative  peak  somewhere  else,  so  that  the  volume 
bounded  by  y  =  0  and  y  =  5*  is  constant.  From  a  hydrodynamic  perspective  the 
present  measurements  do  not  suj^rt  using  the  intermittency-weighted  form  for  the 
displacement  thickness.  On  the  other  hand,  the  assumption  that  the  displacement 
thickness  may  be  modeled  in  this  way  is  valid  for  describing  the  highest-order  acoustic 
source  mechanism.  The  dominant  dipole  source  is  associated  with  the  positive  mass 
flux  deficit  peak,  which  is  a^roximated  by  the  intermittency-weighted  form. 
However,  the  quadrupole  source  is  not  modeled  under  this  assumption. 

5.3  NOISE  FROM  A  NATURAL  BOUNDARY  LAYER  TRANSITION 

This  study  bases  its  conclusions  on  the  behavior  of  single  artificially  generated 
turbulent  spots,  both  in  insolation  and  in  tandem.  In  a  naturally  occurring  transition 
region  the  flow  is  made  up  of  spots  which  at  random  at  the  upstream  end  of 
the  transition  zone,  grow  as  they  convect  downstream,  and  eventually  merge  to  form 
a  low  Reynolds  number  turbulent  boundary  layer.  The  effects  of  spot  intmetion  cm 
spot  growth  have  only  been  partially  considered,  and  the  nature  of  spot  formation  has 
beat  ignored  entirely. 
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It  appears  that  the  degree  of  interaction  rq)iesented  by  the  experiments 
described  here  does  not  modify  the  convection  speed,  the  rise  time,  or  the  mass  flux 
deficit  fluctuation  peak  amplitudes  significantly.  This  result  suggests  that  these 
differences  do  not  modify  substantially  the  order  of  magnitude  estimates  of  the  sound 
level  from  an  isolated  turbulent  spot.  In  addition,  Zilbermann  et  a/.  (1977)  showed 
that  a  turbulent  spot  retains  its  idratity  well  past  the  location  where  it  enters  a 
tuibulrat  environment.  The  data  on  interacting  spots  presented  here  shows  that  for 
a  weak  interaction,  the  poteitial  flow  fluctuations  outside  the  spot  are  not  notably 
atteiuated  from  those  observed  in  the  isolated  ^t  case.  If  the  statements  "the  spot 
retains  its  identity"  and  "the  poteitial  flow  fluctuations  are  not  seriously  modified" 
can  be  seen  as  equivaloit,  thtm  the  effect  of  spot  interaction  on  the  level  of  sound 
radiated  from  the  mature  regions  of  the  transition  zone  is  small  indeed.  The  isolated 
turbulent  spot  measurements,  interpreted  in  the  light  of  the  acoustic  source  model, 
show  that  the  fluid  motion  responsible  for  the  dipole  source  corresponds  to  the  portion 
of  the  turbulent  spot  which  protrudes  above  the  surrounding  laminar  boundary  layo*. 
The  structure  which  Zilberman  et  al.  (1977)  showed  to  persist  well  into  a  tripped 
turbulent  boundary  layer  corresponds  exactly  to  this  structure.  It  may  be  that  the 
early  stages  of  the  fully  turbulent  layer  generate  sound  according  to  the  mechanism 
described  here,  until  the  point  where  the  large-scale  structure  of  the  outer  flow  which 
is  due  to  the  persistence  of  turbulent  ^ts  has  decayed.  Further  exp^moitation 
should  clarify  the  extent  to  which  this  statement  is  correct. 
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To  estimate  the  sound  radiation  from  a  natural  transition  zone  from  the 
radiation  characteristics  of  a  single  turbulent  spot  requires  1)  the  radiation  from  a  spot 
as  a  function  of  location  in  the  transition  zone  and  2)  the  local  spot  density,  i.e. ,  how 
many  qmts  there  are  pm*  unit  area,  again  as  a  function  of  location.  The  extrapolation 
of  the  isolated  spot  results  to  natural  transition  makes  use  of  the  following 
assumptions: 

1.  The  single  spot  radiation  characteristics  are  well-modeled  by  the 
isolated  spot  estimate  given  above.  The  turbulent  spot  radiation,  then, 
retains  its  dipole  character  throughout  the  transition  zone.  This 
assumption  alone  indicates  that  the  transition  r^ion  should  be  a  more 
efficient  sound  radiator  than  the  fully  turbulent  boundary  layer,  which 
radiates  as  a  quadrupole  (Powell,  1960;  Hardin,  1991). 

2.  The  longitudinal  coordinate  in  the  isolated  spot  result  may  be 
replaced  by  ^=(x-X(,)/Ax,  where  Xo  is  the  location  of  the  upstream  edge 
of  the  transition  zone  and  Ax  is  the  transition  zone  length  given 

by  Equation  4.7.  Note  that  this  results  in  the  isolated  spot  radiation 
being  given  for  a  range  of  0.26  <('<  1.30,  extending  beyond  the  end  of 
the  hypothetical  transition  zone. 

3.  The  local  acoustic  source  strength  is  the  local  single  spot  radiation 
weighted  by  the  local  ^t  density,  i.e.,  lNr(x)~p4(x)  N(x)  >  where  p^fx)  is 
the  local  dipole  strength  per  spot  and  N(x)  -N/U  is  the  local  number 
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of  spots  per  unit  length.  N(x)  is  the  local  spot  rate  found  by  Farabee,  et 
al.  (1974)  and  used  by  Lauchle  (1981)  and  Audet  et  (1989a,  1989b)  in 
their  calculations: 

N(f)  =  1.272  (iL)fe-<^"“^’  (5.«) 

Ax 

As  shown  in  Figure  S.4,  the  spot  rate,  N,  rises  from  zero  at  the  upstream  end  of  the 
transition  zone  to  a  maximum  at  about  30%  of  the  transition  zone  length  (see  Lauchle, 
1980),  and  then  back  to  zero  again.  The  total  radiated  pressure  at  any  instant  is  the 
integral  of  p4(x)  k(x)  over  the  transition  zone.  Weighting  the  radiated  pressure  from 
a  single  spot  shown  in  Figure  5.2  with  this  distribution  will  result  in  a  maximum 
located  near  of  the  midpoint  of  the  transition  zone,  as  shown  in  Figure  5.4.  The 
model  of  the  acoustic  contribution  to  the  wall-pressure  used  by  Audet  et  a/.  (1989a, 
1989b)  did  not  consider  the  streamwise  variation  in  AMFD/pU  and  predicted  the 
location  of  the  maximum  acoustic  contribution  to  be  located  at  =  0.5.  The 
streamwise  distribution  of  local  radiated  sound  pressure  is  also  shown  in  Figure  5.4. 
Here,  the  estimate  is  based  on  the  scaling  analysis  presented  above,  with  p,,(x)  given 
by  the  curve  shown  in  Figure  5.2. 

Qualitatively,  the  shape  of  the  streamwise  variation  of  the  acoustic  pressure 
fluctuation  given  by  Audet  et  a/.  (1989a,  b)  is  confirmed  by  wdghting  the  single  spot 
radiation  result  with  the  numb^  of  spots  po*  unit  length.  This  statement  may  be  made 
more  precise  when  the  effect  of  spot  interaction  upon  the  streamwise  variation  in 
radiated  sound  pressure  is  better  understood.  In  any  case,  this  argument  based  on  the 
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Figure  5.4  Extrqralation  of  single  spot  radiation  estimate  to  natural  transition 
ztHie.  Local  acoustic  source  strength  is  prqxnrtional  to  the  product  of  the  local  spot 
source  stroigth,  weighted  by  the  local  spot  density. 


single  spot  behavior  does  indicate  that  the  streamwise  variation  in  the  rms  acousHc 
pressure  shown  by  Audet  et  o/.(1989)  is  largely  due  to  the  spatial  distribution  of  the 
spot  rate.  However,  this  conclusion  is  tentative  at  this  time  because  the  assumptions 
which  have  led  to  it  are  still  largely  unproven,  especially  the  assumption  that  the  local 
single  spot  radiation  remains  dipole-like  throughout  the  transition  zone. 


CHAPTER  6 


CONCLUSIONS 

Measurements  of  the  velocity  field  of  a  single,  isolated  tuifoulent  spot  were 

made  in  order  to  obtain  an  estimate  of  the  ensemble-averaged,  unsteady  local 

fluctuation  in  the  mass  flux  deficit,  and,  in  turn,  the  fluctuation  in  the  mass  flux 

normal  to  the  wall.  These  fluctuatimis  have  been  characterized  by  the  times  for  the 

mass  flux  deficit  to  rise  to  its  maximum,  and  then  fall  to  its  minimum.  In  addition, 

the  streamwise  variations  in  these  fiuctuatimis  were  also  deduced.  The  measured  mass 

flux  deficit  was  also  compared  to  the  form  assumed  in  previous  analyses  of  transition 

noise.  The  componrat  of  velocity  normal  to  the  wall  in  the  freestream,  above  the 

maximum  turbulent  spot  height,  was  also  measured  to  understand  the  streamwise  and 

spanwise  developmoit  of  this  quantity.  Similar  measurements  of  turfiuloit  spot  pairs 

were  carried  out  to  investigate  how  their  interaction  modifies  the  observed  isolated 

spot  results.  The  behavior  of  the  normal  velocity  above  the  turbulrat  spot  was  used 

to  amstruct  a  simple  source  acoustic  model  for  the  turbulent  spot.  An  estimate  of 

the  sound  radiation  from  an  isolated  turbulent  spot  was  obtained  using  a  scaling 

analysis  employing  the  simple  source  model  and  the  scales  of  motion  deduced  from 

€ 

the  expmimental  results. 

6.1  DISCUSSION 

The  mass  flux  deficit  was  calculated  from  measured  ensemble-averaged, 
unsteady  velocity  profiles.  The  shape  of  the  unsteady  mass  flux  deficit  time  series 
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(see  Figure  4.3)  may  be  explained  as  follows.  The  positive  peak  occurs  at  the  same 
location  as  the  maximum  height  of  the  turbulent  spot  which  is  greater  than  the  height 
of  the  surrounding  laminar  boundary  layer.  At  this  location  the  averaged  velocity 
profile  is  identical  to  a  turbulent  boundary  layer.  This  structure  causes  the  part  of  the 
fiow  which  is  initially  above  the  laminar  boundary  layer  to  decelerate,  then  accelerate, 
as  the  turbulent  spot  maximum  height  passes  by,  after  which  the  height  of  the 
turbulent  spot  decreases,  causing  a  decrease  in  the  mass  flux  deficit.  The  minimum 
of  the  mass  flux  deficit  corresponds  to  the  trailing  edge  of  the  turbulrat  portion  of  the 
turbulent  spot.  At  this  point  the  hdght  of  the  turbulait  spot  is  equal  to  that  of  the 
surrounding  laminar  boundary  layer.  Because  the  velocity  profile  is  "fuller"  than  the 
undisturbed  laminar  profile,  the  mass  flux  deficit  is  actually  less  than  that  of  the 
undisturbed  flow.  The  integral  of  the  mass  flux  deficit  over  the  area  of  the  wall  at 
any  instant  is  zero,  indicating  that  the  net  mass  flux  normal  to  the  wall  is  also  zero 
at  any  instant.  The  convection  speeds  of  the  positive  and  negative  peaks  of  the  mass 
flux  deficit  fluctuation  were  in  good  agreement  with  other  investigators’  measurements 
of  the  convection  speeds  of  the  spot  maximum  spot  height  and  trailing  edge, 
respectively.  The  convection  speed  of  the  positive  maximum,  76%  of  the  freestream, 
also  corresponds  well  with  other  measurements  of  the  overall  turbulrat  spot 
convection  speed.  This  quantity  serves  as  a  velocity  scale  for  the  flow  and  is  useful 
for  the  acoustic  analysis.  A  length  scale  for  the  flow  was  obtained  from  the  positive 
peak  amplitude  of  the  mass  flux  deficit. 
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The  time  scale  characterizing  the  change  between  the  fully  laminar  and  fully 
turbulent  flow  was  estimated  by  the  difference  either  in  the  times  of  arrival  of  the 
leading  edge  and  positive  mass  flux  deficit  peak  (the  rise  time),  or  in  the  times  of 
arrival  of  the  positive  peak  and  the  n^ative  peak  (the  fall  time).  Comparison  of  the 
nondimensionalized  estimates  of  this  quantity  at  one  velocity  (40.6  cm/s)  to  a  single 
estimate  obtained  from  one  velocity  profile  at  a  lower  velocity  (19  cm/s)  suggests  that 
the  rise  time  is  indq)endent  of  Reynolds  number,  while  the  fall  time  is  not.  This 
conclusion  finds  support  in  the  results  of  Wygnanski  et  al.  (1982),  which  states  that 
the  leading  edge  and  maximum  spot  hdght  convection  speeds  are  indq)0ident  of 
Reynolds  number,  while  the  trailing  edge  convection  speed  increases  as  Reynolds 
number  decreases. 

The  spatial  distribution  of  the  velocity  fluctuation  nomud  to  the  wall  in  the 
potential  flow,  found  by  direct  measurem^it,  shows  that  the  velocity  peaks  are  in 
phase  across  the  width  of  the  tuibuloit  spot,  and  that  as  the  spot  grows  each  peak 
increases  in  length,  width,  and  hdght,  so  that  for  each  peak  the  mass  flux  normal  to 
the  wall  increases  with  time. 

The  mass  flux  through  a  surface  parallel  to  the  wall,  above  the  maximum 
turbulent  spot  height,  serves  as  the  boundary  condition  for  the  acoustic  wave  equation, 
according  to  the  Liq>mann  Analogy.  Considering  the  mass  flux  through  each  peak 
to  be  an  unsteady  mass  source  or  sink,  the  turbuloit  spot  may  be  described  as  a 
distribution  of  simple  sources.  Because  of  the  low  Mach  number  of  the  flow,  this 
distribution  of  sources  may  be  considoed  compact  since  the  size  of  a  turbulent  ^t 
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is  much  smaller  than  the  wavelengths  of  the  sound  produced  by  the  flow.  The  relative 
magnitudes  of  the  mass  flux  for  each  of  the  peaks  resulted  in  a  simple  source 
distribution  which  could  be  described  as  the  superposition  of  a  longitudinal  quadrupole 
and  a  dipole.  The  axes  of  both  these  sources  are  aligned  with  the  freestream.  The 
dipole  is  associated  with  the  large,  positive  mass  flux  deficit  peak.  The  solution  to 
the  acoustic  equation  given  by  Liq)mann  (1954)  was  approximated  by  a  multipole 
expansion  so  that  the  order  of  magnitude  of  the  contribution  of  each  int^ral  in  the 
expansion  to  the  total  radiated  sound  could  be  determined.  The  variation  in  the  source 
strength  with  streamwise  distance  was  determined  by  substituting  the  scales  found 
experimentally  into  the  multipole  expansion.  These  results  show  that  the  dipole 
contribution  dominates.  The  dipole  stioigth  per  unit  area  scales  as  (AMFD/pU)/^^ 
while  the  dipole  strength  scales  as  AMFD/pU.  The  acoustic  source  is  more  efficient 
in  the  early  stages  of  the  turbulent  spot’s  developmrat  because  of  the  1/t;^  behavior, 
but  the  total  radiated  sound  power  is  less  because  the  size  of  the  turbulent  spot  is 
small  at  this  point  in  its  development.  The  dipole  strength  per  unit  area  decreases  as 
the  spot  grows,  but  the  total  radiated  sound  power  increases  due  to  the  compoisating 
effect  of  the  increased  size  of  the  source. 

In  previous  analysis  of  the  sound  genmted  by  boundary  layer  transition, 
Lauchle  (1981)  and  Lagier  and  Somette  (1984)  assumed  that  flow  quantities  could  be 
described  by  an  intermittency-weighted  form.  Lauchle’s  treatmrat  of  the  displacement 
thickness  was  evaluated  using  the  measured  behavior  of  the  mass  flux  deficit.  It  was 
found  that  the  positive  peak  was  adequately  represented  by  this  model  but  the  negative 
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peak  was  not,  so  that  the  assumed  form  could  not  be  said  to  model  the  hydrodynamic 
fluctuations  accurately.  However,  the  dominant  acoustic  behavior  was  associated  with 
the  motion  which  gave  rise  to  the  positive  mass  flux  deficit  peak,  so  that  the 
intermittency-weighting  assumption  could  be  sera  to  rq>resrat  the  part  of  the  flow 
which  contributes  the  dominant  part  of  the  sound  radiation. 

The  effect  of  turbulent  spot  interaction  on  these  results  was  studied  by 
measuring  velocity  profiles  of  two  turbulrat  spots  generated  at  the  same  location  with 
a  fixed  time  delay,  so  that  the  second  spot  leading  edge  passes  through  the  calming 
region  of  first  spot  and  overtakes  the  trailing  edge  of  the  first  spot.  This  interaction 
scenario  was  studied  briefly  by  Gutmark  and  Blackwelder  (1987),  who  showed  that 
for  a  sufficiently  small  time  delay,  the  convection  velocity  of  the  second  spot  leading 
edge  was  reduced.  The  growth  of  the  leading  edge  of  the  second  turbulrat  spot, 
which  occurs  by  destabilization  of  the  surrounding  laminar  flow,  is  attenuated  because 
the  velocity  profile  is  fuller,  and  therefore  more  stable  in  the  calming  region.  While 
the  present  study  shows  a  significant  modification  of  the  ov^iang  regicm,  the 
convection  velocities  of  the  second  turbulrat  spot  wrae  not  modified  to  any  significant 
degree.  Furthermore,  the  motions  associated  with  the  mass  flux  deficit  positive  peak 
were  also  not  modified.  This  suggests  that  up  to  the  level  of  interaction  investigated 
the  turbulrat  spot  may  still  be  modeled  as  a  dipole  source.  It  may  be  concluded  that 
the  intermittrat  region  of  the  boundary  layer  graerates  sound  more  efficiently  than  the 
fully  turbulrat  r^on  when  the  sqjaration  between  turbulrat  spots  is  large  enough  that 
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the  large  positive  mass  flux  deceit  peak  at  the  ^t  maximum  height  still  exists.  The 
level  of  interaction  for  which  this  behavior  persists  has  yet  to  be  established. 

An  extension  of  the  isolated  turbulent  spot  results  to  a  natural  transition  zone 
was  carried  out  using  two  simple  assumptions:  (1)  the  effect  of  turbulrat  spot 
interaction  on  the  large-scale  bdiavior  of  individual  spots  is  negligible,  and  (2)  the 
local  dipole  source  strength  is  proportional  to  the  product  of  the  isolated  turbulent  spot 
source  strength  and  the  local  turbuloit  spot  density.  Wh^  the  source  strength  was 
scaled  in  this  way  it  exhibited  qualitatively  the  same  behavior  as  that  shown  by  the 
model  of  Audet  et  al.  (1989),  peaking  near  the  midpoint  of  the  transition  zone. 

SUMMARY 

The  conclusions  of  this  study  are  as  follows: 

1.  The  range  of  the  nondimensionalized  rise  time  for  the  isolated  turbulent  spot 
was  found  to  be  0.06<T<0.30  while  the  nondimensionalized  fall  time  spanned 
a  range  0.13<T<0.74.  These  values  confirm  the  lower  range  of  limits  set 

on  this  parameter  by  Lauchle  (1981). 

2.  The  time  scales  describing  the  local  rate  of  change  of  the  streamwise  mass 
flux  deficit  increase  with  distance  from  the  generation  point. 

3.  The  streamwise  spatial  growth  rate  of  the  peak  amplitudes  of  the  mass  flux 
deficit  decreases  with  distance  from  the  generation  point. 

4.  The  peak  amplitudes  of  the  normal  component  of  velocity  increases  with 
distance  from  the  generation  point,  but  levels  off. 


5.  The  results  from  the  isolated  spot  case  were  not  significantly  modified  by 
weak  interaction  between  two  turbuloit  spots  in  the  case  observed. 

6.  The  acoustic  source  due  to  the  large-scale  fluctuation  in  the  turbulent  spot 
flowfield  may  be  modeled  as  the  superposition  of  a  dipole  associated  with  the 
positive  mass  flux  deficit  peak  and  a  longitudinal  quadrupole. 

7.  The  form  of  the  unsteady  displacemoit  thickness  assumed  by  Lauchle  (1981) 
is  invalid  from  a  hydrodynamic  point  of  view  but  adequately  models  the 
dominant  dipole  acoustic  source. 

8.  Extrapolation  of  the  isolated  spot  sound  estimate  to  a  natural  transition  flow 
suggests  that  the  sound  source  strmgth  reaches  its  maximum  near  the  midpoint 
of  the  transition  zone  and  that  the  attire  transition  zone  radiates  as  a  dipole 
source.  It  may  be  concluded  from  the  last  assertion  that  the  natural  transition 
zone  is  a  more  efficient  radiation  source  than  the  fully  turbuloit  boundary 
layer,  which  radiates  as  a  quadrupole. 

6.3  SUGGESTIONS  FOR  FURTHER  WORK 

The  present  study  was  motivated  by  the  need  to  explain  how  the  transition 
region  of  a  boundary  layer  generates  sound  more  efficioitly  than  does  the  fully 
tuibuloit  r^on.  The  focus  of  the  present  work  was  the  bdiavior  of  an  isolated 
turbulent  spot,  the  cause  of  the  intermittent  nature  of  transitional  flow.  The  results 
presented  are  for  an  isolated  turbulent  spot  and  for  a  weakly  intoacting  turbulent  spot 
pair.  Ensuring  the  generality  of  these  results  requires  additional  information 
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concerning  the  effects  of  further  merging  on  individual  turbulent  spot  structure  and 
how  dependent  the  flow  scales  deduced  from  these  measurements  are  on  the  Reynolds 
number.  In  addition,  the  present  results  are  in  some  sense  limited  because  they  are 
based  on  many  realizations  of  point  measurements  taken  over  a  long  period  of  time. 
This  suggests  recourse  to  more  advanced  fiill-field  measurement  techniques.  Finally, 
there  is  a  need  for  direct  measuremrat  of  the  sound  produced  by  both  isolated 
turbulent  spots  and  a  natural  transition  flow. 

The  dominant  sound  production  mechanism  in  a  turbulent  spot  is  the  dipole 
associated  with  the  mass  flux  deficit  positive  peak.  The  extent  to  which  turbulent  spot 
interaction  changes  this  behavior  from  that  observed  in  an  isolated  turbulent  spot  is 
important  for  generalizing  the  single  turbulent  spot  sound  source  model  to  a  natural 
transition  flow.  If  the  qualitative  behavior  of  the  positive  mass  flux  deficit  peak  is  not 
significantly  altered,  then  it  aqrpears  that  the  dominant  sound  source  in  a  natural 
transition  flow  may  indeed  be  seen  as  a  distribution  of  dipoles.  Experimrats  to  test 
the  extent  to  which  this  is  true  should  begin  by  rq)eating  the  interacting  turbulent  spot 
measurements  presented  above  for  shorter  time  delays  between  the  goieration  of  each 
spot.  The  extent  to  which  the  identity  of  a  turbulent  spot  may  remain  distinct  after 
ratering  a  fully  turbulent  environment,  shown  by  Zilbraman  et  al.  (1977),  should  also 
be  examined  with  the  aim  of  detecting  differences  in  the  unsteady  mass  flux  deficit  as 
the  turbulent  spot  passes  by.  These  measurements  would  be  greatly  improved  if  a 
single  realization  of  the  velocity  profile  could  be  obtained  simultaneously,  rather  than 
one  point  at  a  time,  as  in  the  present  measurements.  This  requirement  should  lead  to 
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the  use  of  either  rakes  of  hot  wires  or  hot  films,  or  the  use  of  whole-field  imaging 
techniques  such  as  hydrogen  bubble  visualization  or  Particle  Image  Velocimetry.  The 
apparent  weak  dqiradence  of  the  nondimensionalized  rise  time  and  the  apparent  strong 
dependence  of  the  nondimensionalized  fall  time  on  the  Reynolds  number  should  be 
investigated  more  thoroughly.  This  may  be  done  by  measuring  velocity  profiles  as 
was  done  in  the  present  study  at  differnit  freestream  velocities,  or  with  a  different 
generation  point  location. 

An  important  contribution  to  the  study  of  transition  noise  would  be  the  direct 
measurement  of  the  sound  generated  by  the  intermittent  flow.  This  study  could  start 
with  measurements  of  acoustic  intensity  conditioned  on  the  detection  of  the  presence 
of  an  isolated  turbulent  spot.  These  experiments,  carried  out  over  a  range  of 
freestream  velocities,  would  provide  a  test  of  the  results  of  the  scaling  analysis 
presented  here.  Once  the  clear  relationship  between  the  unit  of  intermittency  (a  single 
turbulent  spot)  and  the  radiated  sound  has  bear  established,  then  measurements  of 
interacting  spots  and  a  naturally  occurring  transition  zone  should  be  carried  out.  A 
study  along  these  lines  is  currently  being  carried  out  in  an  anechoic  facility  at  the 
Applied  Research  Laboratory. 
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